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ABSTRACT 

We present a kinematically-unbiased search to identify young, nearby low-mass mem- 
bers of kinematic moving groups (MGs). Objects with both rotation periods shorter 
than 5 days in the SuperWASP All-Sky Survey and X-ray counterparts in the ROSAT 
All-Sky Survey were chosen to create a catalog of several thousand rapidly-rotating, 
X-ray active FGK stars. These objects are expected to be either young single stars 
or tidally-locked spectroscopic binaries. We obtained optical spectra for a sub-sample 
of 146 stars to determine their ages and kinematics, and in some cases repeat radial 
velocity (RV) measurements were used to identify binarity. Twenty-six stars are found 
to have lithium abundances consistent with an age of ^ 200 Myr, and show no ev¬ 
idence for binarity and in most cases measurements of Ha and v sin i support their 
youthful status. Based on their youth, their radial velocities and estimates of their 
3-dimensional kinematics, we find 11 objects that may be members of known MGs, 8 
that do not appear associated with any young MG and a further 7 that are close to the 
kinematics of the recently proposed “Octans-Near” MG, and which may be the first 
members of this MG found in the northern hemisphere. The initial search mechanism 
was ~ 18 per cent efhcient at identifying likely-single stars younger than 200 Myr, of 
which 80 per cent were early-K spectral types. 
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1 INTRODUCTION 

In a series of publications in the 1960s, Olin Eggen hypoth¬ 
esised the existence of a “Local Associ ation” of co-moving , 
early-type stars in the solar vicinity (lEggenl flOGlI . Il965ll . 
This association (also known as the Pleiades moving group) 
included members of the Pleiades, a Persei and IC 2602 open 
cluster s, and rough ly one-third of the B7-A0 stars within 
300pc llEggenlll98^L 

Following Eggen’s hypothesis, detailed observations 
of several young, chromospherically active, fast-rotating, 
late-type stars led to the suggestion that they too 
were part o f the the Local Association (e.g. AB 
Dor, PZ Tel , Innis. Coates fc Thompso^ 198^ BO Mic, 


lAnders et al.l Il993l : LO Peg. Ijeffries et al 


I994I I . Follow¬ 


up spectroscopy of optical counterparts to coronally ac¬ 
tive late-type stars found in EUV and X-ray surveys 
found that a large fraction were c o-mov ing with the 
Local Association lljeffries fc Jewell Il993ll : many were 
as young or younger than the Pleiades (~ 100 Myr) 


E-mail: a.s.binks@keele.ac.uk 


as evidenced by the presence of lithium in their pho¬ 
tospheres, which is otherwise burned rapidly during 
the pre main sequence in low mass stars l|Jeffriesl 
1995h. Subs e quent work (e.g. Zuckermanfc_Webt 


^rreset^ Aod: Zuckerman et al.l I2001I: iMontes et all 


200 ll : Zuckerman fc Song||2004l l suggested that the Local As¬ 


sociation has kinematic and spatial substructure, consisting 
of several different co-eval, co-moving streams of stars at 
a range of young ages, that have been collectively termed 
“near by young moving groups” (hereafter refer red to as 
MGs. [^ckerrnan fc Songif2004l : [Torres et al.ll2008ll . The ori¬ 
gins of these MGs is still an open question. They may be 
the result of the dissolution of young open clusters, with 
members of each MG sharing a com mon birthplace as wel l 
as common kinematics. For instance, IPe Silva et al.l (l2013ll 
found that stars in the Argus MG and the IC 2391 open 
cluster shared common chemical abundances and kinemat¬ 
ics. However, an abundance analysis of a p opulation of ob- 
jects in the AB Doradus MG (ABDMG) bv iBarenfeld et ^ 
ll2013ll indicates that approximately half of the previously 
suggested members do not share a similar chemical compo¬ 
sition, arguing against a common origin. 
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Finding young stars in the solar neighbourhood is im¬ 
portant because they represent some of the best obser¬ 
vational targets for understanding the early evolution of 
stars and their surroun ding circumstellar environments and 
planetary systems (e.g. IPent et al.ll2013l : iBrandt et al.l[20l3 : 
iBowler et alJuMa l. They are much closer than their equiv¬ 
alents in young clusters and star forming regions, offering 
advantages both in terms of sensitivity and spatial reso¬ 
lution. If stars can be linked to particular coeval MGs, 
then their ages can reasonably be assumed similar to that 
MG as a whole. At ages of 10-100 Myr gas giant planets 
around MG members are expected to be much more lumi¬ 
nous than in older systems, and young stars are frequently 
surrounded by debris discs that may evidence the formation 
of terrestrial planets or provide diagnostic indicators of un¬ 
seen planets. Lower mass MG members potentially provide 
even better targets to investigate planets and circumstel¬ 
lar environments because their lower luminosities enhance 
the contrast with giant planets of a given mass. Examples 
of work that exploits the youth and proximity of MG mem¬ 
bers includes the high contrast infrared imaging detection of 
multiple planets surrounding the 30 My r old A-type dwarf , 
HR 8799, a member of the Golumba MG jMarois et al.llioo^ . 
l2010tl . and the identification of a planet around the AO star 
0 Pic; the eponymous member of the 0 Pictoris MG (BPMG, 
iLagrange et al.lf2010l l. 

Much work has focused on finding new low-mass MG 
members using both kinematic selection and kinematically 
unbiased surveys. Kinematic selection may be efficient at 
discovering new members of known MGs, but precludes the 
discovery of nearby, young objects that are not members of 
these groups. Examples include the proper-motion selected 
searches reported by Schlieder, Lepine & Simon (2010, 2012) 
or the work of Malo et al. (2013, 2014) and Gagne et al. 
(2014, 2015) who used both positions and Galactic veloci¬ 
ties to assign probabilities of membership to new candidate 
members of several known MGs. 

Kinematically unbiased searches are possible but less 
efficient. For example, from a n initial sample of 405 la te-type 
stars within 25 pc of the Sun. lMaldonado et al.l (l2010l i found 
only 6 per cent that may be candidate members of known 
MGs based on their space motions. A more focused approach 
is to pre-select stars which are likely to be young based on 
their magnetic activity. Young stars are magnetically active 
as a result of their fast rotation, convective envelopes and 
consequent dynamo-generated magnetic fields; this activity 
is manifested as chromospheric and coronal emission that 
can be detected via optical emission lines or UV and X- 
ray flux. Exa mples of this a ppro ach can be found in the 
earlier works of Jeffries! lll995h and iMontes et al. (l200lll. bu t 


more recently [Torres et ^ 


Lepine fc Sinionl ] 2009!) . 


d^_Silva_et_a^ ( 20()9h and IShkolnik et al.l ll201l| j. Notably. 


Shkolnik et al.l (l2012l l pre-selected a sample of nearby, X- 


ray active M-dwarfs, finding many new MG members but 
also finding that about 50 per cent of the young M-dwarfs 
could not be assigned to any of the currently known MGs. 

Here we describe a new, kinematically unbiased method 
to select young stars that relies on the fact that stellar ro¬ 
tation is strongly age dependent. At young ages, a large 
fraction of low-mass stars have fast rota tion rates (rota¬ 
tion periods less than a few da ys - e.g. jPatten fc SimorJ 
1 19961 : iKrishnamurthi et al.lIl998h . Angular momentum loss 


due to magnetised stellar winds leads to spin down on a 
mass-dependent timescale, ranging from ~ 50 Myr for G- 
stars to hun d reds of Myr f o r M-d warfs (see for example 
iBarnesllioO^ . lMessina et al.1 (l201(]|j report rotation periods 
of the order of several days or less for many MG mem¬ 
bers, confirming that a selection based on a short rota¬ 
tion period is likely to favour young stars, although may 
be contaminated by members of older, tidally-locked short 
period binary systems. Confirming the youth of candidates 
requires high resolution spectroscopy. Young F- ,G- and K- 
stars should ha ve large abundances of li t hium in their pho¬ 
tospheres (e.g. ISestito fc Randichl l2005l : Ijeffriesi l2006l l and 
multiple radial velocity measurements can be used to iden¬ 
tify short-period binary systems. 

In this paper we report our initial efforts to find young, 
nearby stars using a parent sample of fast-rotating, active 
stars selected from the union of the SuperWASP transit¬ 
ing planet survey and the ROSAT X-ray all-sky survey. 
These candidates were followed up with high resolution spec¬ 
troscopy at the Nordic Optical Telescope (NOT) and Isaac 
Newton Telescope (INT). In 321 we describe how rotation 
periods were estimated from the SuperWASP photometry 
database. The initial candidate selection is described in 33] 
along with details of the spectroscopic observations. In 331 
we describe the techniques used to measure radial and rota¬ 
tional velocities, temperatures, chromospheric activity and 
Li abundances. The multiple methods that were used to con¬ 
strain the ages of the observed targets are discussed in 331 
and in 33] the space motions of the young, Li-rich targets 
are calculated and compared to known MGs. We provide a 
discussion of individual objects in and in 331 we discuss 
the potential of a repeat survey, focusing on the efficiency 
and relative success of this work at identifying kinematic 
sub-structure in the young sample. 


2 PERIOD DETERMINATION 

The SuperWASP project (|Pollacco et al.ir2nn6h is a wide- 
held photometric survey for transiting exoplanets that has 
been operating since 2004. The wide held of view of its two 
instruments located in La Palma (Spain) and South Africa 
(one set of eight Ganon 200 mm f-1.8 lenses on each site, cov¬ 
ering 482 square degrees for each observatory and backed by 
high quality 2048 x 2048 i-Kon GCD detectors) and typical 
observing cadence of ~ 10 minutes also make it prohcient at 
identifying many types of stellar variability with timescales 
from one hour to several weeks. Data are available for tens of 
millions of objects with brightness in the approximate range 
8 < K < 15 covering most of the sky. An initial catalog 
of objects was generated by cross-correlating the ROSAT 
sky-survey (IRXS ) and pointed phase (2RXP) catalogues 
JVoges et al.lll999l . see 1)5.31 for an analysis of the ROSAT 
data acquired in this work) with objects in the SuperWASP 
archive. We created a sample of 5477 stars using the criteria 
that SuperWASP targets must be within either the 3a posi¬ 
tion uncertainty or 10” of the X-ray source (whichever was 
larger), have declinations > —20° and contain more than 
1000 photometric data points in the archive up to the date 
20th July 2010 when the database query was performed. 

Periodic variable stars were identified us ing the Lomb- 
Scargle periodogram technique described in iMaxted et al.l 
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searches for significant periodicities in unevenly sampled 
data. The SuperWASP light-curves are measured over sev¬ 
eral seasons of observation, each of which typically have 
~ 8000 unevenly sampled data points. To measure the pe¬ 
riod, we calculate the normalised power Pn{u}) at a given 
angular frequency, cj = 2ttu. The highest peaks in the cal¬ 
culated power spectrum correspond to candidate periodic¬ 
ities in the time series data. To obtain a solution for the 
light-curve, a least-squares fit of the sinusoidal function 
yi = asin(a;ti) -I- bcos{u}ti) to magnitudes i = 1,2, ...,N is 
found. A power spectrum is obtained based on a chi-squared 
fit of the light-curve: 


where 


PnG) 




Xo 


( 1 ) 


and 







{mi - Vif 


To ensure the signal is not a noise artefact, the false 
alarm probabilities (FAPs) were calculated by making 100 
permutations of each light-c urve using a bootstrap Mont e 
Carlo technique developed bv ICollier Cameron et all ll2009l l. 
The FAP related to a given power P„ is taken as the frac¬ 
tion of randomised light-curves that have a highest power 
peak that exceeds Pn, which is the probability that a peak 
of a given height is merely caused by statistical variations, 
i.e. white noise. The spectrum of FAPs was used to esti¬ 
mate the power value in the periodogram for which this 
probability is 0.1, 1 and 10 per cent. Uncertainties are split 
into 2 separate error bars: first is the (averaged) measure¬ 
ment error using AP = (equation 2 in iMessina et al.l 

I2OIOII . where h/r is the finite frequency resolution of the 
power spectrum and is equal to the full-width half max¬ 
imum of the highest power peak in the frequency spec¬ 
trum. The second error bar (where appropriate) is the 
standard error in 2 or more period measurements. Effects 
such as differential rotation between seasons may exceed 
any error bar generated from the width of the peak fre¬ 
quency, therefore errors b ased on only 1 season are likely 
to be underestimates (e . g. IPeinhold. Reiners fc Basrill2018l : 


lEpstein fc Pinsonneaultl 12014 1. The error bars of objects 

with 2 seasons of data also risk under-estimation because 
of the broad probability distribution of the standard error 
in 2 measurements. 

Based on all the available periodograms, the periods 
were given a quality classification. If 3 or more seasons of 
data for an object resulted in an average period with a stan¬ 
dard error smaller than 5 per cent they were graded as ‘A’. 
To assess the period determination reliability, we measured 
the difference in yf (A^^) between the second highest peak 
in each periodogram to the highest peak. If A^^, averaged 


over all observing seasons, was large enough (each light- 
curve contains ~8000 datapoints, so an average A^^ > 1000 
was chosen as our minimum cut-off) and the value of the 
highest peak was larger than the FAP value at 0.1 per cent, 
we considered this a reliable period measurement. Objects 
with 1 or 2 seasons of data and satisfying the average Pi.yf 
criteria, and had a standard error smaller than 5 per cent 
were graded as ‘B’, whereas objects failing either the /X-yf 
condition or resulting in large standard errors for their pe¬ 
riods were graded ‘C’. 

From the initial catalog of 5477 targets a sub-sample 
of 146 were chosen for spectroscopic investigation based on 
inspection of their light-curves, period analysis, and their 
visibility during the observing run. Rapidly rotating, late- 
type stars with periods less than 5 days were considered 
for target selection. A sub-sample of 26 observed targets 
are listed in Table [T] (all other objects are listed in an on¬ 
line supplement). Note that for the purposes of this paper, 
the light-curves were reanalysed using a larger dataset (see 
and nine of the redetermined periods were subsequently 
found to have periods > 5 days. 

We provide the raw light-curves, the periodograms and 
corresponding phase-folded light-curves for all objects that 
we obtained spectroscopy for (see 1)3.l|l as an online supple¬ 
ment to this paper. A summary of the targets investigated 
in this paper and their measured periods are provided in 
Table [ 1 ] along with their uncertainties and quality classifica¬ 
tion. Note that the light-curves and the periods quoted in 
these tables are based on a reanalysis of all the SuperWASP 
data available up until December 2012, performed after our 
spectroscopic follow-up. In most cases these agree with the 
original periods which were based on data up to the 20*^ 
July 2010, however there are a handful of significant dis¬ 
crepancies. 

The main source of B and V magnitudes are 
from the AAVSO All-Sky Photometric Survey (APASS, 
iHenden et al.ll2012ll . which covers the magnitude range 10 < 
V < 17. All objects with V magnitudes < 10 that were 
unavailable in the APASS catalog were instead sourced 
from the NOMAD catalog (lZacharias_etalJ 200^. All K 
magnitudes are from 2MASS ( Cutri et al.l l20o'^ Proper- 
motions are extracted fro m the PPMXL catalog of posi¬ 
tions and proper motions (IRoeser. Demleitner fc Schilbachl 
I2OIOII . BVK photometry and proper-motions are presented 
in Tables |4] and [ 6 l respectively. Figure [T] displays the entire 
observed sample plotted in terms of their rotation periods 
and V — Ks colour. 


3 SAMPLE SELECTION AND 

SPECTROSCOPIC OBSERVATIONS 

3.1 Sample selection 

Prior to selecting targets for spectroscopy, the light-curves 
for the sample available in the SuperWASP catalog were 
checked to identify objects most likely to be showing vari¬ 
ability in their light-curves as a result of starspot modula¬ 
tion. These were more likely to be objects with several sea¬ 
sonal light-curves. Whilst some objects had consistent peri¬ 
ods measured from many seasons of data, with clear sinu¬ 
soidal behaviour, some of the sample had sparse, poor qual¬ 
ity data and sometimes inconsistent period determinations. 
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Figure 1. Rotation periods and V — Ks colour for the entire observed sample. Red squares represent 26 objects that were later identified 
to be likely-young and were considered likely-single stars (referred to as the ‘likely-young sample’, see EDI- All these objects are labelled 
from A to Z and their corresponding ROSAT IRXS names are provided in Table ^ Any objects that had two or more separate RV 
measurements that varied by more than 5kms~^ were assumed to be close-in short-period binaries and are denoted here by open green 
circles. All other objects are designated filled green diamonds. This symbol sche me for all other plots will remain the same throughout this 
work. The various lines represent the I- and C-sequence gyrochrones proposed bv lBarn^ l|2Q03h and calibrated bv iMamaiek fc Hillenbrandl 
ll2008lj . The method of inferring an age range from the gyrochrones is described in ani- note the three objects A, B and L which are 
used as examples to describe how to fit the age ranges. 


Because of the range of quality in data the selection pro¬ 
cess was entirely subjective. Any objects indicating a sharp 
dip in their light-curve over any part of their phase were 
flagged as eclipsing binaries and discarded. The light-curve 
variability amplitudes were restricted to < 0.2 mag and ob¬ 
jects with measured periods to within a thousandth of ei¬ 
ther a half or integer number of days were discarded because 
they risked incorrect period determinations due to aliasing 
effects. Priority was given to objects where the profile of the 
multi-seasonal light-curves significantly changed in shape; a 
strong indication of starspot migration over the course of 
several months. 

The range of target spectral-types is restricted to the 
range mid-F to early M (approximately 1.0 < V — Kb < 3.5) 
for two reasons. First, because lithium is used as a primary 
youth indicator (see section ESI, this range contains stars 
where significant and measurable Li depletion is e xpected 
on timescales of 10 — 100 Myr (see ISoderblomll2010l and ref¬ 
erences therein). The hot boundary is where Li-depletion 
timescales become very long, and poorly understood diffu¬ 
sive or oth er non-standard mixing pro cesses may become 
dominant (iBoesgaard fc Tripiccd Il986li . The cool bound¬ 
ary marks the point where Li depletion occurs extremely 
rapidly and even ve ry young stars may have depleted their 
Li fe.g. I Jeffrie jl200(^ L Secondly, low-mass MG stars will usu¬ 
ally rotate very rapidly; FGK stars typically have spin-down 
timescales of ~ 50 — 100 Myr, and M-dwarfs ~ 300 Myr or 


longer (this is illustrated by the gyrochrones in Figure [T]). 
Thus fast rotation in M-dwarfs is not necessarily a sign of 
youth and the targetting of such objects could lead to high 
contamination with older stars. The basis for choosing our 
objects in terms of colour and rotation period is discussed in 
more detail in HS.ll where we use ‘gyrochronology’ to provide 
crude age estimates. 

Objects in the SuperWASP catalog may have correlated 
with ROSAT sources by random chance. To test this, the 
original positions of all objects observed in the ROSAT cat¬ 
alog were offset by 2’ and checked again for neighbouring 
SuperWASP sources. The ROSAT error circles are about 
10 — 20”, therefore a search radius of 15” was used to iden¬ 
tify any neighbouring objects in the SuperWASP catalog. 
Only one object (target W) out of 146 had a neighbouring 
SuperWASP source subsequent to being shifted by 2’. The 
chance of random correlation is 0.68 per cent. 

An important point is that the nature of the target 
selection, from an all-sky survey, allows for a kinematically 
unbiased sample, independent of proper-motion or radial ve¬ 
locity (RV) criteria. This provides an opportunity to identify 
young stars that do not share the galactic space motions of 
previously identified MGs. The sample is not spatially un¬ 
biased, because the SuperWASP survey avoids galactic lat¬ 
itudes between ±20°. 


© 2015 RAS, MNRAS 000,[TH22] 































Searching for New Nearby Young Moving Groups in the Northern Hemisphere 5 


Table 1. Measured periods for the sample of 26 likely-young objects (see SMli. Periods are the reanalysed values (and standard e rrors) 
using the analysis in [[2l The first error bar is the (averaged) uncertainty calculated using equation 2 in iMessina et aP (l2010ll . The 
second error bar (where appropriate) is the standard error in measurements for 2 or more seasons. ^N’ refers to the number of seasonal 
lightcurves analysed for each target star. The column labelled ‘Q’ refers to the quality of the period determination. The ‘Date’ column 
refers to the month and year when the first observation of the target in this work took place, Jll = June 2011, and D12 = December 
2012. X-ray count rates (CR) and HRl ratios are extracted from either the IRXS or 2RXP catalogs and logLx/Lhol is the X-ray to 
bolometric luminosity ratio. The equivalent data for all other objects observed in this work (and for all subsequent tabular data) are 
available online. 


SuperWASP ID/Label 
(ISWASP J-) 

V 

(mag) 

V-K, 

(mag) 

Period 

(days) 


N 

Q 

Date 

CR 

HRl 

log Lx/+bol 

012457.96+255702.4, 

A 

10.718 

1.811 

3.048+0.120+0.026 

1711 

2 

B 

D12 

0.06 

0.42 

-3.35 

013514.32+211622.4, 

B 

10.788 

1.490 

1.871+0.010+0.004 

1994 

3 

A 

D12 

0.01 

1.00 

-3.84 

023503.81+313922.1, 

C 

10.337 

1.743 

1.274+0.013+0.002 

2038 

4 

A 

D12 

0.15 

0.17 

-3.18 

030405.14+300309.6, 

D 

10.875 

1.807 

1.806+0.020+0.005 

1518 

2 

B 

D12 

0.06 

0.66 

-3.23 

031628.14+563857.7, 

E 

10.667 

2.103 

3.850+0.080 

584 

2 

C 

D12 

0.08 

0.34 

-3.31 

032231.55+285319.8, 

F 

10.672 

2.013 

1.654+0.016+0.011 

281 

2 

C 

D12 

0.06 

0.86 

-3.34 

034319.02+222657.2, 

G 

11.571 

2.545 

2.893+0.036 

1075 

1 

B 

D12 

0.06 

1.00 

-3.07 

050206.19+311102.2, 

H 

10.792 

2.218 

4.531+0.116 

4169 

1 

B 

D12 

0.07 

-0.07 

-3.48 

052146.83+240044.4, 

I 

10.393 

2.323 

3.436+0.054 

9955 

1 

C 

D12 

0.30 

-0.08 

-3.01 

084748.63+342356.8, 

J 

10.304 

1.686 

4.231+0.074 

996 

2 

C 

D12 

0.03 

-0.08 

-3.90 

133708.20+444454.5, 

K 

11.050 

1.954 

4.376+0.196 

301 

1 

C 

D12 

0.07 

-0.11 

-3.32 

135458.37-054354.1, 

L 

9.231 

2.091 

5.361+0.140+0.040 

319 

2 

C 

D12 

0.09 

-0.60 

-4.16 

143854.57+330019.9, 

M 

11.036 

2.219 

5.095+0.156+0.036 

3576 

2 

B 

D12 

0.05 

0.40 

-3.36 

162506.55+300225.8, 

N 

10.205 

3.022 

1.002+0.006+0.098 

51786 

2 

B 

D12 

0.22 

0.29 

-3.27 

214537.36+271110.8, 

O 

11.425 

2.305 

1.267+0.015+0.225 

1152 

3 

B 

D12 

0.06 

0.29 

-3.23 

133241.69+223006.6, 

P 

9.655 

2.808 

1.040+0.013 

186 

1 

C 

Jll 

0.84 

-0.13 

-2.97 

155007.38-022211.5, 

Q 

10.050 

2.538 

2.241+0.019+0.007 

37160 

3 

A 

Jll 

0.48 

-0.09 

-2.99 

162946.59+281038.0, 

R 

10.654 

1.925 

1.426+0.011+0.003 

5023 

4 

A 

Jll 

0.09 

0.87 

-3.14 

171808.56+250612.0, 

S 

10.614 

3.027 

2.420+0.047+0.008 

5432 

2 

B 

Jll 

0.27 

0.04 

-3.08 

172228.64+365842.1, 

T 

10.504 

2.800 

1.229+0.011+0.001 

10864 

9 

A 

Jll 

0.25 

-0.06 

-3.14 

173103.32+281506.5, 

U 

10.213 

2.402 

1.263+0.012+0.002 

13828 

3 

A 

Jll 

0.29 

-0.05 

-3.11 

180426.56+393047.1, 

V 

11.532 

2.365 

1.547+0.017+0.004 

3260 

5 

A 

Jll 

0.07 

0.01 

-3.15 

192502.00+442950.7, 

W 

9.939 

2.398 

1.205+0.013+0.001 

11141 

2 

B 

Jll 

0.24 

0.18 

-3.23 

205830.74-090223.3, 

X 

11.180 

2.071 

2.251+0.026+0.012 

5173 

6 

A 

Jll 

0.11 

0.28 

-2.99 

225617.59+205236.2, 

Y 

11.410 

2.731 

1.100+0.009+0.003 

1164 

4 

A 

Jll 

0.08 

0.06 

-3.21 

230752.70+171015.2, 

Z 

10.780 

2.881 

1.073+0.007+0.002 

980 

3 

C 

Jll 

0.30 

-0.25 

-3.01 


3.2 Observation strategy 

We obtained high-resolution multi-echelle spectroscopy of 
146 objects over 2 telescope runs using the Fibre-fed Echelle 
Spectrograph (FIES) on the 2.56 m Nordic Optical Telescope 
(NOT) at the Roque de los Muchachos Observatory. Our 
first observing run of 68 targets lasted 4 consecutive nights 
from 21 June 2011, and a further 78 targets were observed 
on 27, 28 and 29 December 2012. FIES was used in medium 
resolution mode (A/AA = 46,000). Each target observation 
was bracketed with short ThAr arc lamp exposures to pro¬ 
vide a wavelength calibration. In what follows, objects ob¬ 
served in June 2011 or December 2012 are denoted Jll and 
D12, respectively (these are listed in the ‘Date’ column in 
Table [1} • 

Our observing strategy was to obtain a single observa¬ 
tion of a target, reduce the data in real time at the telescope 
and inspect the spectrum for the presence of the Li i 6708A 
doublet. If Li was clearly present then, if necessary, a fur¬ 
ther exposure was performed to obtain a signal to noise ratio 
(SNR) per pixel of ^ 50 around the Li line. Objects where 
Li was detected were observed again on a subsequent night 
in order to check for short-term RV changes that might be¬ 
tray their close binary nature (see 44.5l for further details on 
binary contaminants in the survey). 

The real time data reduction at the telescope and sub¬ 


sequent reduction at Keele University were accomplished us¬ 
ing the FIEStool package created by Eric Stempels. This 
package performed bias, blaze and flat field corrections, op¬ 
timal extraction of the spectra and wavelength calibration 
using the adjacent arc lamp exposures. With fibre bundle 
#3, FIES covers a spectral range of AA 3630-7260A over 78 
spectral orders. Observations were also made of several RV 
and low-activity template stars (selected for their extremely 
low levels of chromospheric Call H and K emission) over a 
similar spectral-type range as the targets. These were used 
to calibrate RVs and projected rotational velocities (wsinz) 
and are listed in Table [2] 

To further identify whether or not objects were short- 
period binaries, several of the most probable young objects 
were later re-examined using long-slit spectroscopy on the 
Isaac Newton Telescope (INT) on the 23'^'^ and 24*’*' March 
2013, in combination with the H1800V grating and IDS- 
235 wide-held camera. This allowed for an RV precision of 
~ 1 — 2kms“* and these measurements were compared to 
the RVs measured at the NOT. In this work, the RV mea¬ 
surements at the INT are not incorporated into the hnal 
averaged RV, rather they are used to check for RV consis¬ 
tency and Li content. All individually measured RVs are 
presented in column 3 of Table |3] and INT observations are 
subscripted with an ‘i’. 
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Table 2. RV and min im um a ctivity standards - All RV 
are from IChu bak et alj ||2012|) . except for HD 114762 
jUdrv. Mayor ^.^Quel^ Il999h ~ All V sin i values are from 


Glebocki &: Gnaciriski| |2005h . All logi?jjj^ are from 


Isaacson Sz Fisched ll20] 

except 

for HD 

3651 JPacd l2013h 

and HD 190007 jMittae 

et alJl2011^ 



Standard 

SpT 

Ry ^Rv 

kms“^ kms“^ 

nsini 
km s“^ 

log Rhk 

HD 114762 

F8V 

49.4 

0.5 

1.7 

-4.902 

HD 1461 

GOV 

-10.158 

0.116 

5.0 

-5.008 

HD 95128 

GIV 

11.293 

0.108 

3.1 

-4.973 

HD 217014 

G2.5IVa 

-33.118 

0.128 

2.8 

-5.054 

HD 197076 

G5V 

-35.402 

0.030 

2.9 

-4.872 

HD 9407 

GOV 

-33.313 

0.124 

0.0 

-4.986 

HD 115617 

G7V 

-7.844 

0.128 

0.5 

-4.962 

HD 101501 

G8V 

-5.464 

0.105 

2.3 

-4.483 

HD 3651 

KOV 

-32.940 

0.042 

0.6 

-4.849 

HD 10780 

KOV 

2.814 

0.086 

0.9 

-4.700 

HD 131509 

KOV 

-44.749 

0.167 

6.4 

-5.142 

HD 4628 

K2.5V 

-10.229 

0.030 

1.6 

-4.979 

HD 190007 

K4Vk 

-30.270 

0.105 

2.8 

-4.592 

HD 209290 

MOV 

18.275 

0.120 

3.8 


HD 119850 

M1.5V 

15.778 

0.061 

1.8 


GJ 411 

M2.0V 

-84.683 

0.030 

1.6 


GJ 526 

M4.0V 

15.778 

0.061 

1.8 



4 DATA ANALYSIS AND RESULTS 
4.1 Equivalent widths (EWs) 

Li EWs were estimated from the NOT spectra using the 
SPLOT procedure in IRAF0 Lines were measured by man¬ 
ual integration, where the linear continuum is subtracted 
and pixels were summed with partial pixels at the ends 
to obtain a flux. Although the Li line at 6707. SA actually 
consists of 2 neighbouring lines at AA 6707.761, 6707.912A 
our spectroscopy cannot resolve these lines. There is also a 
contaminating weak Fel spectral line at 6707.441 A blended 
into the Li resonance line. Where necessary, we use a cor¬ 
rection for t his blend, Wa( b le nd) = 20(B — V) — 3mA 
provided in ISoderblom et all lll993al l. EW measurement 
errors were estimated using the procedure highlighted in 
ICavrel de Strobel fc Snitj lll988l 'l: 

JEW = 1.6 X \/FWHM X p/SNR, (2) 

where FWHM refers to the observed full width half max¬ 
imum (in A) of the line, and p is the pixel size (0.113A). 
The SNR was calculated empirically from the rms of fits to 
continuum regions. 

Measurements of the 6563A Hq line EW were made by 
direct integration above a surrrounding pseudo-continuum. 
We provide both Li and Her EW measurements in Table |4] 


^ IRAF is distributed by the National Optical Astronomy Ob¬ 
servatory, which is operated by the Association of Universities 
for Research in Astronomy (AURA) under cooperative agreement 
with the National Science Foundation. 


4.2 Radial velocities 

We used the FXCOR procedure in iraf, which cross¬ 
correlates pixels between the target object s and those of 
temp late spectra to determine a relative RV llTonrv fc David 
Il979h . SNR improves towards the red end of the spectrum, 
therefore measurements were made over 9 consecutive or¬ 
ders (AA 5920-6520A), avoiding the broad Na D lines and 
contamination from telluric absorption. To determine the 
centre of the CCF peak a Gaussian fit was made to pixels 
within the top 40 per cent of the normalised peak. Some or¬ 
ders match better than others and poor CCF matches could 
be due to a spectral-type mismatch between an RV stan¬ 
dard and a target, too much noise in the data, or, if the 
target is a close binary, then the CCF may contain multiple 
(unresolved) peaks. The spectrum of a binary companion 
could be blended and more apparent in some orders than 
others. In some cases, no CCF peak was observed at all. In 
order to provide a more robust RV measurement (and also 
assess the precision) the clipped mean RV from 9 orders. 
Any measurements > 2a from the mean were removed, and 
the process iterated a maximum of ten times. A weighted 
average relative RV for each target was calculated over each 
order for each t empla te star using the R quality factor in 
iTonrv fc David (Il979l 1^. In all cases the template that pro¬ 
vided the lowest error in RV was used. Finally the RVs were 
transferred onto a heliocentric reference frame. 

There were three sources of error: the standard error 
measured over the echellogram orders (weighted by R), the 
published RV uncertainty of the standard star used in the 
calibration and a systematic uncertainty estimated from the 
cross-correlations of all the RV standards with one another. 
These errors were treated as independent and were added 
in quadrature. RV measurements and errors are provided in 
Table 131 along with any previous literature values. The RV 
measurements for each RV standard relative to one another 
formed a matrix of cross-correlation values. The systematic 
uncertainty was measured to be 0.36 kms“^ for the tem¬ 
plates used in the June 2011 run and 0.27 kms”*^ for the tem¬ 
plates used in the December 2012 run. To avoid any spectral 
mismatch between target and template calibrations, tem¬ 
plates are restricted to within 0.5 spectral classes of the tar¬ 
get star in each case. RVs obtained from INT spectroscopy 
were measured using a similar procedure, but this time sep¬ 
arating the long-slit spectra into ten 50A windows between 
6600 and 7100A and obtaining the mean RV, weighted by 
R. These measurements were not included in the final av¬ 
eraged RV however because the overall accuracy was much 
poorer, limited by systematic uncertainties of ~ lkms“^ in 
the wavelength calibration, were used only to support the 
data obtained from the NOT. 

From the initial sample of 146 targets, 28 were found to 
have indeterminate RVs (as a result of poorly constrained 
cross-correlation peaks due to noisy data, high v sin i or per¬ 
haps binarity) and 14 objects were observed to have RVs 
that varied by more than Shms”*^ on the timescale of the 
observing run, presumably as a result of binarity (see 114.511 . 


^ This is the ratio of the height of the true peak in the GCF 
compared to the average CGF value over the whole spectral order. 
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Table 3. Radial velocity measurements. Individual RVs are pro¬ 
vided in column 3 and the final averaged RV in column 4. The 
INT RV measurements are subscripted with an ‘i’ and are not 
factored into the final averaged RV for each target. The column 
titled is the binary score based on the criteria described in 

El 


Target 

HJD 

(2450000-I-) 

-^k^ndiv 

(kms-1) 

(kms-1) 

B 

A 

6289.422 

7.2 ± 0.7 

7.3 ± 0.7 

1 


6290.409 

7.4 ± 0.7 



B 

6289.358 

5.0 ± 1.9 

5.0 ± 2.0 

1 


6290.427 

5.0 ± 2.2 



C 

6289.370 

0.3 ± 2.5 

0.3 ± 2.4 

1 


6290.442 

0.2 ± 2.3 



D 

6289.476 

9.8 ± 2.4 

9.8 ± 2.4 

1 


6290.448 

9.8 ± 2.4 



E 

6289.488 

-2.0 ± 0.6 

-2.1 ± 0.6 

1 


6290.460 

-2.1 ± 0.6 



F 

6289.498 

8.2 ± 1.8 

9.0 ± 1.8 

1 


6290.469 

9.9 ± 1.9 



G^- 

6291.569 

6.0 ± 0.5 

6.0 ± 0.5 

1 

H 

6291.589 

11.2 ± 0.6 

11.2 ± 0.6 

1 


6376.436 

14.2 ± 2.2; 



Y 

6291.598 

13.5 ± 0.5 

13.5 ± 0.5 

1 


6376.429 

17.1 ± 2.2; 



J 

6289.580 

10.1 ± 0.6 

10.2 ± 0.6 

1 


6290.617 

10.2 ± 0.6 




6375.474 

11.6 ± 1.5; 



K 

6289.720 

-6.4 ± 0.6 

-6.4 ± 0.6 

1 


6290.708 

-6.4 ± 0.6 




6375.619 

-1.2 ± 3.5; 



L 

6289.731 

-4.1 ± 0.6 

-4.2 ± 0.6 

1 


6290.719 

-4.2 ± 0.6 




6375.631 

-2.3 ± 1.6; 



M 

6289.738 

-19.2 ± 0.6 

-19.2 ± 0.6 

1 


6291.634 

-19.3 ± 0.6 



N 

6289.780 

-68.0 ± 0.7 

-68.5 ± 0.7 

1 


6290.791 

-69.1 ± 0.7 



O 

6290.317 

-22.5 ± 1.8 

-22.9 ± 1.6 

1 


6291.314 

-23.4 ± 1.4 



P 

5735.502 

-5.1 ± 0.7 

-5.1 ± 0.7 

1 


6375.594 

-5.0 ± 2.2; 



Q" 

5736.457 

-10.3 ± 0.3 

-10.6 ± 0.4 

1 


5737.503 

-10.9 ± 0.4 




6375.688 

-10.4 ± 2.3; 



R 

5735.567 

-25.3 ± 0.8 

-25.2 ± 0.8 

1 


5736.515 

-25.1 ± 0.7 




6375.699 

-20.3 ± 2.2; 



S 

5736.483 

-34.7 ± 0.3 

-34.8 ± 0.3 

1 


5737.513 

-34.9 ± 0.3 



T 

5735.613 

-24.6 ± 0.6 

-25.3 ± 0.6 

1 


5736.524 

-26.0 ± 0.6 



U 

5735.643 

-18.2 ± 0.5 

-18.2 ± 0.5 

2 

V 

5736.591 

-25.5 ± 0.6 

-26.4 ± 0.6 

1 


5737.569 

-27.3 ± 0.7 




5734.615 

-33.5 ± 0.6 

-33.5 ± 0.6 

1 

X 

5736.652 

-3.0 ± 0.4 

-2.7 ± 0.4 

1 


5737.617 

-2.4 ± 0.4 



Y 

5737.676 

-2.5 ± 0.6 

-2.5 ± 0.6 

2 

Z 

5737.662 

-8.4 ± 0.3 

-8.4 ± 0.3 

2 


a: RV = 6.06 zb 0.29kms“^ jMermi lli od. Mayor &: Udrvll2009h . 
b: RV = 14.2 di 1.4k ms~^ ijBia zzo et~ al.l l2012ll , c: RV 
= —13.6 zb 1.8kms~^ jKordopatis et " | 2013ll , d: RV = 
—33.1, —33.1, —32.2 km s~^ llFrasc^^ alJ |201lD~ 


4.3 Projected rotational velocities 

Projected rotational velocities (rsini) were estimated from 
the FWHM of the CCF using a set of simulations for objects 
observed on the NOT. Broadened versions of a number of 
slowly rotating template stars (v sin i templates, see Tabled 
were produced by convolving high SNR spectra (~ 100) with 
a rotational broadening convolution kernel, using a limb- 
darkening parameter e = 0.6. Because the usinf templates 
all had a small, but finite, rotational speed (< 5kms“^) a 
correction was made for the small measured v sin i values 
of the standards to set their non-rotating profiles at zero, 
and artificial Gaussian noise was injected into the spun-up 
spectra. For stars with large rotational broadening, the high- 
frequency (wavenumber k > 400 m”*^) components in the 
spectra were filtered out to minimise the effects of Gaus¬ 
sian noise (and the same filtering also applied to the target 
spectra). Each standard was broadened in 1 kms“^ steps be¬ 
tween 5 and TOkms”*^ and the FWHM of each broadened 
standard cross-correlation peak was recorded using FXCOR 
and repeated for the same echellogram orders used in the 
RV analysis. Relationships between FWHM of the GCF and 
V sin i were derived using the means from the orders weighted 
by the R quality factor (see lTonrv fc Davislll979lf for each 
template. To ensure there were no spectral mismatches, only 
usinf template spectra within 5 spectral sub-classes of the 
target were used in each case. The resolving power and SNR 
for spectra on the INT were insufficient to obtain v sin i mea¬ 
surements. The measured v sin i values are presented in Ta¬ 
ble IH] 

4.4 Effective temperatures and Li abundance 
calculations 

Effective temperatures were estimated based on a fifth 
order polynomial interpolation of the V — Ka and Tefr 
values between FO and M 3 stars provided in table 5 of 
IPecaut fc MamaiekI (l2013lj and spectral-types are calculated 
from their linear interpolation with V — Ka in the same 
table. There is evidence to suggest that the metallicities 
of you ng, nearby stars do not vary much from the solar 
value llSpina et al.ll201^ and FGK stars should be close to 
the ZAMS by 20 Myr. The EW of the Li l 6708A line is 
transformed into logarithmic Li abundances (on the usual 
scale where logV(H)= 1 2) using the curves of growth from 
ISoderblom et all lll993bl f and correcting for non-local ther- 
modynamic eq u ilibriu m effects using the code provided in 
ICarlsson et al.f (Il994h . Table |4] displays the temperatures, 
EWs, Li abundances and their associated uncertainties. 


4.5 Binary contaminants 

Tidally-locked, short-period binaries (herein, TLSPBs) 
maintain high rotational velocities well into their main se¬ 
quence lifetimes and therefore can retain short rotation peri¬ 
ods into old age. According to numerical calculations, solar- 
mass binaries starting from the birthline that have periods 
< 8 days will have circular orbits an d synchronous rotatio n 
upon reaching the main sequence (IStahler fc Pallal l2005l f . 
They will not spin down with age like single stars, but 
their rapid rotation will ensure they maintain their mag¬ 
netic activity levels. Furthermore there is evidence that TL- 
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Table 4. EWs, BVK photometry, temperatures, spectral-types (SpT) and Li abundances. Li EWc is the final EW after correcting for 
the blended Fe I line. 


Target 

Li EW 

Li EWc 

Ha EW 



B 



V 



Ks 

Teff 

SpT 

A(Li) 


(mA) 


(m 

A) 

(A) 


(mag) 


(mag) 


(mag) 

(K) 



A 

223 

213 

± 

13 

1.28 

11.470 

± 

0.045 

10.718 

± 0.042 

8.907 

± 

0.019 

5430 

G9 

2-89+^"« 

B 

154 

144 

± 

21 

2.22 

11.387 

± 

0.008 

10.788 

± 0.020 

9.298 

± 

0.016 

5850 

G1 

2.9©°;1“ 

C 

209 

197 

± 

20 

0.99 

10.980 

± 

0.061 

10.337 

± 0.049 

8.594 

± 

0.017 

5510 

G 8 

9 

^•^^-0.15 

D 

193 

186 

± 

21 

1.77 

11.579 

± 

0.087 

10.875 

± 0.076 

9.068 

± 

0.018 

5440 

G 8 

9 77+0.14 
' * -0.18 

E 

222 

211 

± 

13 

0.82 

11.482 

± 

0.049 

10.667 

± 0.023 

8.564 

± 

0.020 

5120 

K1 

^•^^- 0.12 

F 

184 

174 

± 

25 

0.71 

11.356 

± 

0.050 

10.672 

± 0.059 

8.659 

± 

0.019 

5210 

K1 

9 Ko+O.lS 
^•^^-0.18 

G 

267 

246 

± 

25 

0.49 

12.537 

± 

0.056 


11.571 

9.026 

± 

0.020 

4710 

K4 

9 oc+O.lO 
^•'^'^-0.18 

H 

249 

238 

± 

16 

0.80 

11.671 

± 

0.078 

10.792 

± 

0.017 

8.574 

± 

0.021 

5000 

K2 

^•'^-^-0.13 

P 

386 

366 

± 

14 

-0.05 

11.367 

± 

0.045 

10.393 

± 

0.043 

8.070 

± 

0.021 

4900 

K3 

Q 1 q+0.09 
'J--Lo_o.14 

J 

165 

155 

± 

15 

1.83 

10.958 

± 

0.023 

10.304 

± 

0.038 

8.618 

± 

0.022 

5590 

G 6 

'^-0.14 

K 

260 

245 

± 

10 

1.07 

11.852 

± 

0.027 

11.050 

± 

0.022 

9.096 

± 

0.018 

5270 

KO 

^•oy_o.ii 

L 

169 

159 

± 

12 

1.98 

10.089 

± 

0.144 

9.231 

± 

0.104 

7.140 

± 

0.020 

5130 

K1 

9 qo+0.15 

^.oo_o_i9 

M 

106 

93 ± 8 

0.96 

11.934 

± 

0.055 

11.036 

± 

0.048 

8.817 

± 

0.020 

5000 

K2 

-‘-•y^_0.i4 

N 

86 

69 ± 15 

-0.40 

11.311 

± 

0.032 

10.132 

± 

0.037 

7.183 

± 

0.027 

4710 

K5 


0 

254 

242 ± 24 

-0.08 

12.014 

± 

0.032 

11.279 

± 

0.037 

9.120 

± 

0.018 

5060 

K2 

2 . 68 l°;“ 

pb 

159 

142 ± 14 

-0.62 

10.669 

± 

0.043 

9.655 

± 

0.024 

6.847 

± 

0.023 

4390 

K5 

-*-•'^'-0.15 

Q 

174 

155 ± 23 

- 0.01 

10.939 

± 

0.041 

10.050 

± 

0.058 

7.512 

± 

0.023 

4490 

K5 

1 ck+ 0-13 
J-.DO_o.i 8 

R 

208 

196 ± 18 

1.14 

11.394 

± 

0.062 

10.654 

± 

0.041 

8.729 

± 

0.020 

5300 

KO 

'^-0.15 

S 

135 

117 ± 4 

- 1.22 

11.533 

± 

0.184 

10.614 

± 

0.069 

7.587 

± 

0.020 

4330 

K 6 

1 Qi +0.12 

J-.Dl_o.i 2 

T 

62 

44 ± 7 

-0.34 

11.601 

± 

0.046 

10.504 

± 

0.039 

7.704 

± 

0.017 

4500 

K5 

-*-•^-*--0.16 

U 

278 

265 ± 36 

0.19 

11.069 

± 

0.026 

10.213 

± 

0.011 

7.811 

± 

0.016 

4830 

K3 

^.O0_o.25 

V 

269 

252 ± 36 

-0.18 

12.603 

± 

0.066 

11.532 

± 

0.046 

9.167 

± 

0.017 

4860 

K3 

z.oo_Q 26 

W‘= 

222 

208 ± 20 

0.01 

11.011 

± 

0.047 

9.939 

± 

0.061 

7.541 

± 

0.018 

4840 

K3 

9 91 ■l"0.12 
^.D1_o.i7 

X 

337 

279 ± 42 

0.76 

12.208 

± 

0.252 

11.180 

± 

0.112 

9.109 

± 

0.021 

5150 

K1 

9 Qq+0.23 
Z.yo_Q 29 

Y 

141 

122 ± 8 

-0.36 

12.390 

± 

0.071 

11.410 

± 

0.072 

8.679 

± 

0.016 

4560 

K4 

1.D^_0.i7 

Z 

163 

146 ± 30 

-1.50 

11.867 

± 

0.054 

10.780 

± 

0.024 

7.899 

± 

0.027 

4440 

K5 

1.DU q 21 


a: Ho EW = —2.14 A; Li EW = 390 mA l|Li__^Jiu||l99d). Ho EW = —4.00 A; Li EW = 350 mA, b: Ho EW = 0.6 A dMason et ahlllOOhli . 
c: logg = 4.39; [Fe/H] = —0.01 llPinsonneault et al.ll2012f) . 


SPBs may dep lete their lithium at a different rate to sin¬ 
gle s t ars (e.R. iBarrad o v Nava scues. Delivannis fc Sta uffeij 
I2001I : ICanto Martins et al. I2OIII : Strassmeier et al. 2012ll . 
For these reasons we attempt to identify these objects in 
our sample and when we find them we do not attempt to 
estimate their age. 


TLSPBs are expected to be the main contaminants in 
our survey, but should have large RV variations detectable 
by two observations separated by ~ 24 hours. In general, the 
RVs of likely-single, slow rotators with a reasonable SNR 
could be measured to a precision of ~ Ikrns”^. We set a 
criterion that an RV difference > 5kms“^ between 2 mea¬ 
surements of the same object on separate nights identifies 
these objects as close binary systems. A numerical grading 
system is used to distinguish likely single stars (allocated a 
grade of 1) from those very likely to be binaries (allocated 
a grade of 5). A score of 5 was given if we detected RV dif¬ 
ferences > 5kms“^ for a target on separate nights or (if 
there was only one measurement) if there were literature 
sources indicating that either the object is a close binary or 
report a RV measurement > 5kms“^ discrepant from our 
measurement. In addition, to score 5, the average error in 
RV measurement must be less than 5kms“^. Objects are 
graded 4 if there is only one spectrum which results in a) 
a clear, multi-peaked cross-correlation function (CCF) and 
b) an RV error < Skras”*^. A grade of 3 was given if the 
status of the star from the CCF was unclear (presumably 


as a result of poor SNR and/or large usinf), resulting in ei¬ 
ther an indeterminate RV or a RV uncertainty > 5kms“^. 
A grade of 2 was used for objects which had a single spec¬ 
trum, an RV uncertainty < Shms”*^ and a distinct single 
peak in the CCF. Finally, objects scoring 1 had consistent, 
low-uncertainty RV measurements for 2 or more spectra. In¬ 
dividual RV measurements and binary scores for each object 
are presented in Table [S] 

From the entire observed sample of 146 objects, 14 
scored 5, 29 of them scored either 4 or 3, 45 scored 2 and 
29 objects scored 1. Any objects that scored 5 were auto¬ 
matically assumed to be close binaries. All targets selected 
for further analyses in have binary scores of either 1 or 
2. However, 2 or more consistent RV values for an object is 
not a guarantee that the object is not a TLSPB. Given that 
some of the objects in this sample have orbital periods on 
the order of ~ 24 hours, consistent RV measurements could 
be a result of measuring two points at the same phase on an 
RV curve. 


4.6 Binary simulation 

Although the criterion that two separate observations result¬ 
ing in RV differences > 5 kms“^ is sufficient to flag these ob¬ 
jects as TLSPBs, a fraction with RV differences < 5kms“^ 
may remain undetected. In a TLSPB system the rotation pe¬ 
riod will be equal to the orbital period. Given that these ob- 
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servations are incapable of distinguishing between TLSPBs 
and single stars, a simulation was carried out to calculate 
the probability that each object showing no RV variations 
(i.e. not > 5kms“^) was in fact a short period binary seen 
at low inclination or at two similar phases in its orbit. 

The simulation used 10^ iteratio ns, a flat mass-ratio be¬ 
tween 0.20 and 0.95 llHoeevee ^ Il990l) and the primary mass 
was estimated from its spectral-type. Tidally-locked binaries 
will be circular, therefore the eccentricity is zero in all cal¬ 
culations. A value between 0 and 27r was randomly assigned 
for the argument of periastron. A second RV measurement 
was calculated at a later time equal to the cadence (cad) 
between observations. 

We carried out the simulation for all objects that had 
2 separate RV measurements at the NOT. In some cases, 
the separation between the observations was similar to 
an integer number of rotation periods. In these cases, a 
TLSPB nature cannot be ruled out with any confidence 
but for the majority of cases we are ~ 90 per cent con¬ 
fident that they cannot be TLSPBs. In other words in 

> 90 per cent of the simulations there would have been 
a ARV > 5kms“^. Simulations for two objects, target 
P (P = 1.04 days, cad = 2.00 days, the probability that 
the object is a TLSPB (Ptlspb) = 0.12) and target D 
(P = 3.85 days, cad = 0.97 days, Ptlspb = 0.15) are pre¬ 
sented in Figure [2] In the initial target selection, light-curves 
indicative of eclipsing binaries were filtered out, therefore it 
is unlikely that high inclination, eclipsing TLSPBs have been 
included. For example, a 0.7Mq/0.5Mq binary system with 
an orbital period of 3 days would be at least partially eclips¬ 
ing if the inclination angle, i > 60°. Given a random orien¬ 
tation of orbital rotation axes, the probability of observing 
an object with i > 60° is 0.5, therefore Ptlspb could be 
over-estimated by a factor of 2. 

The incidence of close binaries in a random sample 
of field stars is expected to be quite low. We simulated a 
field population based on the binary fraction, period and 
ecc entricity distribution o f field solar-type stars proposed 
by iRaghavan et al.l ll2010ll . finding that only 2.2 per cent 
of such a sample should lead to observed RV differences of 

> 5kms“^ in observations separated by ~ 1 day. However, 
of the 55 objects that we observed twice, 11 were found to 
show evidence of close binarity - significantly higher than 
2.2 per cent. The reason is likely that we are not observ¬ 
ing random field stars, but a sample of fast-rotating, active 
stars, which will tend to b e either young or in close binary 
systems fe.g. IJeffrieslll995l L 

The field star simulation served an additional purpose. 
Where we identify a Li-rich, potentially young star, with 
two consistent RV measurements (within 5kms“^), then as¬ 
suming that these stars share the same binary properties as 
average field stars, there is still a 7.2 per cent probability 
that the true systemic velocity could be different to that 
measured by > 5kms“^ due to the presence of an unseen 
or unresolved wide binary companion. 


5 FINDING YOUNG STARS 
5.1 Gyrochronology 

Gyrochronology provides a colour-dependent method to ob¬ 
tain an age estimate for single, late-type stars by using em- 
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Figure 2. Simulations to calculate the fraction of tidally-locked 
binary objects that would result in an RV difference < 5kms“^ 
over the cadence of two observations for target D (top) and tar¬ 
get P (bottom). The input parameters and a description of the 
simulation is provided in EM The left panel demonstrates the 
probability that a binary star has a particular RV difference of a 
certain value. The right panel is the cumulative probability that 
an object has IRV 2 — RVi| less than a given value. 


pirical fits to rotation rates in open clusters of known age. 
Figure [ 1 ] displays the rotation periods against V — Ks for the 
entire observed sample, overplotted with curves for the in¬ 
terface (I) and convective ( C) sequences of the form d efined 
by in Barnes (2003; 2007). [Mamalek fc Hillenbrandl (l2008fl 
find age discrepancies between the calibrated gyrochrones 
in Barnes (2003; 2007), claiming that the Barnes calibra¬ 
tions would underestimate the Hyades age. We theref o re use 
the updated calibration in iMamaiek fc Hillenbrandl (l2008fl 
an d transform B — V to V — by interpolating table 5 
in IPecaut fc MamaiekI ll2013l) . At this point we stress two 
issues a bout our use of gyrochronology: firstly, the calibra¬ 
tions in IMamaiek fc Hillenbrandl (l2008l ) are for stars with 
0.5 ^ {B — V)o ^ 0.9 (~ F7-K2) and roughly half of our 
sample are redward of this range. Secondly, we find that in 
almost all of our sample the upper age li mits are not g en¬ 
erally affected by the ch o ice of either the iBarM i 1 I 2007 I) or 
IMamaiek fc Hillenbrandl (I 2 OO 8 I) calibrations. 

The bifurcation point is the point at which the G se¬ 
quence meets the I sequence and the isochrones become bi- 
modal redward of this point. To estimate an age using gy- 
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rochronology the following analysis is carried out (see Fig¬ 
ure [T]): 

•If the object lies above the I sequence, it is definitely 
older than this gyrochrone. 

•If the object lie below the I sequence and blueward 
of the the corresponding C sequence, then this target is 
younger than that gyrochrone. 

•If the object lies below the I sequence and redward 
of the corresponding C sequence, there is an ambiguity and 
one must look for an older gyrochrone which satisfies the 
previous condition. 

To clarify how gyrochronology works in practice we 
highlight three stars from Table [T] (targets A, B and L): 
Target A - Target A lies below the 100 Myr I sequence, but 
above the 50 Myr I sequence. It is redward of the 10, 50 and 
100 Myr C gyrochrones, but lies blueward of the 300 Myr 
C gyrochrone and the 300 Myr bifurcation point. Hence the 
age of target A is estimated to be < 300 Myr, with a lower 
limit of 50 Myr. 

Target B - Target B has a minimum age of 50 Myr, as 
it lies above the I sequence of this isochrone. It is below 
the I sequence of all gyrochrones with ages > 100 Myr but 
is only blueward of the C sequence for those with ages ^ 
100 Myr. Hence the age of B is estimated to be between 50 
and 100 Myr. 

Target L - Target L is at least older than 100 Myr, because 
it is above the I sequence at this age. Although it lies below 
all the I sequences ^ 100 Myr it is only blueward of the C 
sequence at an age of 300 Myr. Hence L has a gyrochronology 
age between 100 and 300 Myr. 

The limitations of the initial target selection in terms 
of colour and period now become clearer when examining 
Figure [T] A period < 5 days implies an age < 700 Myr for 
F-stars; < 300 Myr for G-stars with 1.5 < V — Ks < 2.0, but 
then increases again to < 700 Myr for K-stars with 2.0 < 
V — Ks < 3.0 and even older for M-dwarfs. Reversing this 
argument, to be sure an object has an age ^ 100 Myr would 
require it to have a period of < 1 day if it were a late F- 
star, < 2 days for a G-star and then less than 0.7 of a day 
for K-stars and less than 0.3 of a day for M-stars. To avoid 
ambiguity with interpreting the G-sequence we provide (and 
employ) only upper-age limits from gyrochronology for the 
likely-young sample and these are provided in Table (5] 

Gyrochronology will not work for tidally-locked binary 
systems. Because this is the likely status of all the objects 
flagged as 4 or 5 in the binarity tests, we choose not to 
assign ages to them based on gyrochronology (or activity - 
see section ISTl) . 


5.2 Li EW measurements compared to known 
clusters 

The strength of photospheric lithium in cool stellar atmo¬ 
spheres can be used as an empirical age indicator. Stars 
are born with the Li abundance of the interstellar medium, 
but it is depleted by proton capture in their interiors. Mix¬ 
ing processes bring the Li-depleted material to the surface, 
resulting in an age- and mass-dependent photospher ic Li- 
abundance (see lSestito fc Randichl[2005l : ljeffriesll2006ll . 

Whilst theoretical models partially capture the be¬ 
haviour of PMS Li depletion in FGK stars, they are strongly 
sensitive to the assumed opacities in the atmospheres and 


interiors of these stars and cannot readily explain the spread 
in Li abundance that is seen at a given effective tempera¬ 
ture in presumably coeval clusters. For example, the spread 
in Li abundance amongs t K-stars in the Pleiades is 1 — 2 dex 
dSoderblom et aLlll993d l. For these reasons Li is used here 
only as an empirical age estimator by comparing the Li 
EWs of targets, as a function of B — V and V — Ks (or 
equivalently, Li abundance vers us Teff), with stars obs erved 
in the Hyades (age ~ 625 Myr, [^ rryman et al.||l 998ll and 
Pleiades (age ~ 125 Myr, IStariffer. Schultz fc KirkpatrickI 
Il998ll clusters. Li EW data obtained from Soderblom et 
al. (1990,1993a,b,c,1995a,b); Jones et al. (1996,1997); Jones, 
Fischer & Soderblom (1999); Soderblom et al. 1999 and 
IWilden et al.l (120021 '). In addition, comparisons were made 
with three younger clusters; NGC 2264 (age 5 My r, iKind 


1 19981 : ISoderblom_el^jL _ 1^99 ; Dahm fc Simon! 12005 ), 7 Vel 
(age 10 My r, Jeflhles_et_^ 201-^ and IG 2602 (^ 3 0 Myr, 
iMeola et al.ll2000l : iRandich et al.ll200ll : lRandicbll200ll l . 

In Figure [ 3 ] we present the Li EWs and colours for the 
entire observed sample, and ages are estimated using com¬ 
parisons to the Li patterns in the aforementioned open clus¬ 
ters. One object (labelled as target I) with duplicate, con¬ 
sistent RV measurements has a Li EW/K — Ks indicative 
of stars younger than the Pleiades (more similar to the pat¬ 
tern observed in 7 Vel at ~ 10 Myr). A further twenty-five 
targets with binary scores of 1 or 2 have Li EWs/colours (or 
Li abundances/temperatures) consistent with the Pleiades 
or IG 2602. The remainder of the targets appear to have de¬ 
pleted their Li sufficiently that they appear older than the 
Pleiades, or they are close binaries and their Li-depletion 
may not properly reflect their true age. There is some evi¬ 
dence that close binary system componen ts have anomalous 
Li-de pletion histories (see, for example IStrassmeier et ahl 
l2012li . The bottom plot of Figure [3] is A(Li) versus Tefi. As 
this is essentially the data from the top two plots of Figure [3] 
folded through a transformation, it contains no new infor¬ 
mation on the ages, but does illustrate that the most Li-rich 
stars found are still depleted from an assumed initial level 
of about A(Li)=3.3, so are probably not very young PMS 
stars. The estimated Li age ranges for the entire observed 
sample are presented in Table [5] 


5.3 X-ray activity in the ROSAT catalog 

The initial selected targets described in Sj3]were chosen based 
on their short rotation periods and assessment of apparent 
rotation modulation in their light-curves. All such objects 
originated from a catalog of objects that had entries in both 
ROSAT and SuperWASP. Throughout the selection process 
no criteria were used based on any specific X-ray property 
- a database entry in either the IRXS or 2RXP catalog was 
deemed sufficient. Data were available for all the selected 
targets and were analysed subsequent to the telescope ob¬ 
servations. Lx/Lboi was calculated for each object using the 
following formula: 


Ax ^ _/x_ 

Lboi ^ 2.48 X 10-5 X 


(3) 


The ROSAT catalogs were revisited to obtain count 
rates, exposure times and hardness ratios (GR, exp and 
HRl, respectively, all provided in Table[T]). X-ray fluxes (/x) 
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V-K 



Figure 3. Top: Li EWs and corresponding V — Ks colours. Green 
open circles and filled diamonds represent the same objects de¬ 
scribed in Figure ^ Objects are compared to members in the 
Hyades (625Myr), the Pleiades (125Myr), IC 2602 (30Myr) and 
7 Vel (10 Myr, see text for references) to estimate an Li-based age 
range. Error bars are not included, but are provided in Table [4| 
Bottom: Li abundances as a function of surface temperature. Both 
the target sample and all ancillary data from the open clusters 
are folded through the same curve of growth (described in ^4.4II . 
Fourth-order polynomial fits are overplotted for the Hyades and 
Pleiades, however these are not implied to represent the trend 
of Li abundances in these clusters. We display only the objects 
from Tabled in the bottom panel and Li abundances for all other 
objects are provided in the online supplementary data. 
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Figure 4. log Lx/^bol as a function o{ V — Ks- The 10**' to 90*** 
percentile range of Lx/Tbol values in NGC 2547, the Pleiades 
and the Hyades are represented by the light blue/blue/dark blue 
lines, respectively. 


were calculated by multiplying the cou nt rate by the energy 
conve rsion factor (EOF) provided in IStelzer fc Neuhauseil 
ll200lll : 

ECF = (8.31 + 5.30HR1) x lO'^^ergcm"^ cr"\ (4) 

where cr“* is the number of counts per second. Bolo- 
metric magnitudes were found using V magnitudes and 
a main-sequence bolom etric correction inte r polat ed from 
V — Ks using table 5 in IPecaut fc MamaiekI J 2 OI 3 II . Main- 
sequence conversions should be appropriate for the F, G 
and K stars in this sample, however, if stars were younger 
than 30 Myr this may result in a small calibration error of 
~ 0.1 dex. 

Figure 0] shows Fx/Fboi as a function oi V — Ks- For 
reference, the 10**' and 90**' percentile Fx/iboi range for F, 
G and K stars in NG C 2547, the Pleiades and the Hyades 
(based on figure 12 in Ijeffries et al.l I 2 OO 6 I I are represented 
by light blue/blue/dark blue lines (respectively). Fx/Fboi 
is a relatively crude age indicator and largely rotation- 
dependent in any case. Almost all of the candidates could 
be consistent with a Pleiades age (as expected for their ro¬ 
tation rates). As it turned out the probable young objects 
were slightly more active than the average star in our spec¬ 
troscopic sample. 


5.4 Ha emission 

A potentially more direct activity-age dependent probe is to 
observe the strength of the Ha line at 6563A. Ha emission 
lines are diagnostic of strong magnetic activity in the pho¬ 
tospheres of stars, which are linked to rotation and young 
ages. 


3.5 
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The empirical behaviour is that young stars show 
Ha emission, but this emission is not apparent in stars 
with photospheres warmer t han a temperature that ap¬ 
pears to be age-d ependent llReid. Hawlev fc Mated 1 19951 : 
iHawlev et al.lll9^ b Using data from clusters ranging from 
30-625 Myr (IC 2602 IC 2 391, NGC 2516, Pleiades and 
Hyades) iHawlev et al.l (Il999l l derived a log-linear fit to the 
age (in Myr) of a cluster and the V — I colour at which stars 
first display Ha emission given as: 


(U - I). 


emission 


= -6.42 -h 0.99(log age/Myr) (5) 


At redder colours there is a point where no stars 
have Ha in absorption. Compiling data in 7 open clusters 
aged between 30 and 625 Myr (N GC 2391 and NGC 2602, 


SUmffer et al.l 1 19971 : NGC 2547. Ijeff ries. To tten fc JamesI 


2000|; Blanco 1, Panagi fc 0’Del| 1997 ; Pl eiades 


Soderblo m et al.l Il993d~ N GC 2516, IHawlev et abl 1 19991 : 
Hyades. IStaufferet al.l 19911 1 we derived a relationship for a 
given V — I and age at which a star has reduced in magnetic 
activity and starts to show Ha in absorption: 



1.0 1.5 2.0 2.5 3.0 3.5 

V-K, 


{V - /)absorption = -4.19 -|- 0.68(log age/Myr) (6) 

These relationships can be used to estimate the 
minimum/maximum age of a star with an Ha absorp¬ 
tion/emission line. Targets with Ha EW > 200 mA are con¬ 
sidered as emission lines (providing a maximum age) and 
EW < -200 mA for Ha in absorption (a minimum age). 
Targets with Ha EWs in the range ± 200 mA were consid¬ 
ered as ‘filled-in’ lines and equations [5] and |6] were used to 
provide a likely age range. In Table |4] the Ha EWs for the 
likely-young sample are presented, and the ages estimated 
for the entire sample are provided in Ta ble [H V — Ig colours 
are co nverted to V — Ks using table 5 in IPecaut fc MamaiekI 
ll2013ll . The Ha EWs of the entire observed sample are plot¬ 
ted as a function of their V — baseline in Eigure[5] 


5.5 Radii-colour isochrones 

We calculated projected rotational velocities {v sin i) follow¬ 
ing the procedures in i)4.3l and multiplied them by their ro¬ 
tational periods to obtain measurements of the projected 
stellar radius Rsini = 0.02Pnsini (in solar radii). Be¬ 
cause of the sini ambiguity, this technique is only capa¬ 
ble of providing lower limits to the radius. Values of Psini 
are co mpared to radius/V — Ka i sochr ones from the mod¬ 
els of ISiess. Dufour fc ForestiiiH (l2000f l. Between ages of 
10 — 100 Myr a star is contracting on the PMS, and be¬ 
cause the radii calculated from Psini are minima, this sets 
an upper limit to the age. The P sin i values used correspond 
to the Icr lower limit based on both the period and vsini 
measurement, which sets a likely upper age limit. Psini 
offers good age discrimination when dR/dt is large, but be¬ 
comes poor when dR/dt ~ 0 as objects approach the ZAMS. 
The targets that displayed no Li yet appear to have large 
Psini values are more likely to be binaries as opposed to 
being genuinely young. One object in Figure [6l target M, 
has an Psini age (5 — 10 Myr) in contradiction with the Li 
age (30 — 200 Myr) and is at odds with a Ha lower limit of 
32 Myr. 


Figure 5. Ho EW as a function of V — Kb- Objects displaying 
emission may be compared to an upper age limit based on chro¬ 
mospheric activity and a lower age limit for Ho absorption. Only 
objects with EWs > 200 mA (< —200 mA) were identified as true 
emission (or absorption) lines. 



V-K3 


Figure 6. Radii-colour relationship for the entire observed sam¬ 
ple (minus the six objects with an indeterminate usini measure¬ 
ment). Co ntinuous lines represent isochron es generated from the 
models in ISiess. Dufour &; Forestinil i|20QOlj . The plotted points 
are Rsini values and therefore lower limits to the true radius. 
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Table 5. Age estimates in Myr for the entire observed sample, 
based on Gyrochronology, Ha, Li EW and i?sin j/colour. The fi¬ 
nal age estimate is solely from the Li EW age, other age indicators 
are used only as supporting evidence for the Li age. 


Target 

Gyro 

(Myr) 

Li EW 
(Myr) 

Ha 

(Myr) 

Rsini 

(Myr) 

A 

<300 

30-150 

>21 


B 

<100 

30-200 

>17 


C 

<100 

30-200 

>20 


D 

<300 

30-200 

>21 


E 

<300 

30-200 

>27 


F 

<300 

30-200 

>24 


G 

<500 

30-150 

>41 


H 

<500 

30-150 

>32 


I 

<500 

5-30 

35-50 


J 

<300 

30-150 

>19 


K 

<300 

30-150 

>22 


L 

<300 

30-150 

>27 


M 

<500 

100-200 

>32 

<5 

N 

<500 

30-150 

<189 


0 

<300 

30-150 

29-39 


P 

<300 

30-150 

66-127 


Q 

<500 

30-150 

41-64 


R 

<300 

30-150 

>22 


S 

<700 

30-150 

<233 

<80 

T 

<500 

100-200 

<123 


U 

<300 

30-150 

>37 


V 

<300 

30-150 

<53 


w 

<300 

30-150 

37-55 

<80 

X 

<300 

30-150 

>26 


Y 

<300 

30-150 

<99 


Z 

<500 

30-150 

<158 



5.6 Combined age estimation and identification of 
young, single stars 

In terms of providing a secure identification of youth, we 
regard Li as the primary age indicator. Alternative age in¬ 
dicators are used only to support an age determination. If 
no Li is present in the spectra, then the object is automat¬ 
ically assumed to be at least older than the Hyades and is 
assigned a nominal age of 1 Gyr (for the purpose of a paral¬ 
lax estimation - see asm. 

Subsequent analysis in this paper is restricted to 26 ob¬ 
jects that have Li-derived ages younger than 200 Myr and 
have binary scores of 1 or 2. We categorise these objects as 
the likely-young sample. Table [5] shows that, with the excep¬ 
tion of target I, all of the likely-young sample have Li-based 
ages which overlap with the ages constrained from Ha. Each 
age estimate for the likely-young sample is provided in Ta¬ 
ble [5] (all other objects are listed online). Rsini values are 
only identified as younger than main-sequence for 3 objects, 
two of which (targets S and W) are consistent with Li and 
Ha and the other, target M, is discussed in il5.5l Target I 
has an age range of 5-30 Myr, significantly younger than the 
rest of the sample. 


6 KINEMATIC ANALYSIS 
6.1 Parallax estimation 

Only four objects from the entire observed sample have 
a previously measured trigonometric parallax, only one of 


which is part of the likely-young sample. For the rest of the 
objects, parallaxes were estimated using V — Kb (assum¬ 
ing no reddening) a nd a maximum and minimum ag e (see 
Table [5]). Using the ISiess. Dufour fc Forestinil (l200oh evo¬ 
lutionary models (with solar metallicity and no convective 
overshoot) a maximum and minimum absolute K magnitude 
was calculated by interpolating Mk for the youngest and 
oldest isochrone corresponding to the estimated age range, 
from which a ‘photometric’ parallax was calculated. 

In Figure[3the Mk range is plotted for the likely-young 
sample, along with all other objects and TLSPBs. The plot 
shows that the younger objects span a larger range in Mk 
due to PMS contraction. For example, target I, with an 
age range from 5 to 30 Myr has a possible Mk range of 
~ 1.5 mag, whereas target V, with an age ranging between 
30 and 70 Myr results in a difference in Mk of ~ 0.5 mag. 
Redward ot V — Kb = 3.0 isochrones are separated by no 
more than 0.7 mag between 30 Myr and the ZAMS and no 
more than 0.3 mag between 50 Myr and the ZAMS. 

To test the precision of the photometric parallaxes, ab¬ 
solute K magnitudes for the 4 objects with trigonometric 
parallaxes were estimated in the same way; one of which, 
target P, is a member of the likely-young sample, and is 
plotted as a purple downwards-pointing triangle. For each 
measurement the photometric parallax appears to under¬ 
predict the trigonometric parallaxes. The average difference 
between the Mk measured from trigonometric parallax and 
estimated Mk is ~ 0.4 mag and the average range in Mk 
from photometric parallaxes is 0.44 mag. The parallax range 
for all objects is provided in column 6 of Table |51 

6.2 Calculating galactic space velocities 

Using the prescription in Ijohnson fc SoderblomI lll987ll . 
Galactic space velocities (and their errors) were computed in 
terms of U, V and W (where U points towards the Galactic 
centre, V in the direction of Galactic rotation and W in the 
direction of the Galactic north pole). The photometric paral¬ 
lax range is used if no trigonometric parallax is available. Ta¬ 
ble |6] lists all of the required input parameters and columns 
5, 6 and 7 are the calcul ated UVW. Proper motions are 
from the PPMXL catalog (iRoeser. Demleitner fc Schilbachl 
I2OIOII . Two separate error bars are quoted, the first includes 
contributions from the and ctrv uncertainties and 

the second corresponds to half of the range in each veloc¬ 
ity coordinate resulting from the extrema of the possible 
photometric parallaxes (provided in column 6). Should a 
trigonometric parallax exist, or if an object is assigned a 
single-valued main-sequence age (IGyr) then all errors are 
incorporated into one single error bar. 

In Figure [H]the likely-young sample are plotted on Boet- 
tlinger U — V and V — W diagrams. Their space motions are 
compared with 10 nearby MGs (the UVW used to define 
one sigma errors for the MGs are taken from table 1 in 
iGagne et al.ll20i3 and are listed in Table [Til. With the ex¬ 
ception of target P, which has a measured parallax, velocity 
error bars on this plot do not include an error due to the 
parallax uncertainty. Instead, a line connects UVW points 
calculated at the extrema of the distances inferred from the 
photometric parallaxes. The object plotted in blue is target 
I which had an age range significantly younger than the rest 
of the sample (see Table 0. We also note that our analysis 
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Table 6. Kinematic data for the likely-young sample. 


Target 

(mas yr“^) 

(mas yr“^) 

TT 

(mas) 

(kms 

U 

V 

(kms“^) 

W 

(kms“^) 

V sin 2 
(kms“^) 

A 

11.9 ± 1.6 

-14.0 ± 1.7 

7.22-8.31 

-7.4 ± 0.9 ± 0.2 

-4.6 ± 0.9 ± 0.6 

-10.2 ± 0.9 ± 0.4 

8.0 ± 1.1 

Ba 

12.2 ± 1.5 

-0.6 ± 1.5 

5.60-5.62 

- 10.6 

± 

1.5 

-3.9 ± 1.5 

-1.9 ± 1.6 

35.5 ± 0.5 

c 

38.9 ± 1.4 

-21.9 ± 1.4 

8.08-9.04 

-12.9 ± 1.9 

± 

0.7 

-21.1 ± 1.3 ± 1.2 

-1.5 ± 1.3 ± 0.1 

20.5 ± 1.1 

D 

11.9 ± 1.7 

-4.2 ± 1.7 

6.67-7.64 

-12.3 ± 2.1 

± 

0.3 

-3.2 ± 1.4 ± 0.5 

-2.4 ± 1.4 ± 0.1 

18.0 ± 0.5 

E*" 

27.5 ± 2.1 

-35.1 ± 2.1 

9.48-11.75 

1 

o 

b 

H- 

p 

G 

± 

1.2 

-16.2 ± 0.8 ± 1.6 

-6.9 ± 1.0 ± 0.7 

8.6 ± 0.5 

F 

12.4 ± 1.6 

-14.6 ± 1.7 

8.82-10.65 

- 10.2 ± 1.6 

± 

0.2 

-5.8 ± 1.0 ± 0.8 

-5.5 ± 1.0 ± 0.2 

14.6 ± 0.6 

G'= 

22.5 ± 2.1 

-46.2 ± 2.1 

8.78-11.37 

- 6.6 ± 0.6 

± 

0.2 

-21.9 ± 1.0 ± 3.0 

-10.9 ± 0.9 ± 1.1 

6.2 ± 0.5 

H 

27.2 ± 1.3 

-52.5 ± 1.3 

9.76-12.34 

-13.7 ± 0.6 

± 

0.3 

-23.8 ± 0.6 ± 2.9 

-5.6 ± 0.6 ± 0.5 

5.4 ± 0.2 

I"! 

10.0 ± 1.2 

-48.0 ± 1.2 

6.46-12.74 

-11.5 ± 0.5 

± 

0.6 

-25.4 ± 0.7 ± 8.2 

-11.9 ± 0.7 ± 3.3 

6.5 ± 0.4 

J 

-17.1 ± 2.2 

-50.6 ± 2.1 

7.87-8.55 

-12.5 ± 0.9 

± 

0.2 

-29.5 ± 1.2 ± 1.2 

-5.3 ± 1.1 ± 0.5 

5.8 ± 1.2 

K 

-16.9 ± 1.6 

7.4 ± 1.6 

7.06-8.39 

-10.7 ± 1.0 

± 

0.9 

-5.2 ± 1.0 ± 0.3 

-5.5 ± 0.7 

5.4 ± 0.5 

L 

-115.3 ± 1.4 

-69.2 ± 1.5 

18.20-22.50 

-15.1 ± 0.4 

± 

1.4 

-27.6 ± 0.4 ± 3.0 

-5.5 ±0.5 ± 0.2 

4.3 ± 1.2 

M 

-35.7 ± 1.9 

-17.5 ± 1.9 

11.02-11.08 

-8.9 

± 

0.8 

-21.4 ± 0.8 

- 11.0 ± 0.6 

20.4 ± 1.3 

N 

39.4 ± 1.6 

-37.7 ± 1.6 

25.25-33.18 

-25.0 ± 0.4 

± 

1.0 

-37.2 ±0.4 ± 0.1 

-52.6 ±0.5 ± 0.7 

8.1 ± 0.9 

O 

16.6 ± 1.7 

-7.5 ± 1.7 

7.45-9.42 

-8.3 ± 1.0 

± 

0.5 

-23.7 ±1.5 ± 0.3 

-1.3 ± 1.1 ± 1.0 

14.6 ± 0.5 

pe 

-135.0 ± 1.1 

55.2 ± 1.1 

21.71 ± 1.64 

- 20.2 ± 0.2 

± 

0.9 

-5.4 ± 0.2 ± 0.2 

-1.4 ± 0.7 ± 0.2 

19.6 ± 0.4 

Q 

-62.8 ± 1.2 

-72.5 ± 1.2 

17.57-22.84 

-6.9 ± 0.3 

± 

0.2 

-23.7 ± 0.3 ± 3.0 

-5.5 ±0.3 ± 0.1 

16.2 ± 0.5 

R 

28.2 ± 1.4 

-12.4 ± 1.4 

8.25-9.76 

1 

CO 

b 

H- 

p 

± 

0.8 

-8.1 ± 0.7 ± 0.5 

-28.8 ± 0.8 ± 1.0 

22.3 ± 0.4 

S 

-4.0 ± 1.7 

76.3 ± 1.6 

21.69-28.66 

-32.1 ± 0.3 

± 

1.7 

-14.7 ± 0.3 ± 1.0 

-13.4 ± 0.3 ± 0.6 

16.1 ± 0.5 

T 

-37.3 ± 2.5 

63.4 ± 2.4 

24.10-24.19 

-22.2 

± 

0.5 

-18.0 ± 0.5 

-5.7 ± 0.5 

16.8 ± 0.5 

U 

5.0 ± 1.5 

63.2 ± 1.7 

14.71-18.86 

-25.0 ± 0.5 

± 

1.9 

-3.1 ± 0.5 ± 1.2 

-5.4 ±0.5 ± 0.4 

18.3 ± 0.5 

V 

-22.4 ± 3.1 

4.1 ± 3.1 

7.78-9.95 

- 10.8 ± 1.6 

± 

0.1 

-0.6 ± 1.1 ± 27.0 

-0.2 ±1.6 ± 1.4 

18.8 ± 0.5 

w 

30.9 ± 1.5 

7.1 ± 1.5 

16.61-21.26 

-12.1 ± 0.4 

± 

0.6 

-29.2 ± 0.5 ± 0.3 

-13.8 ± 0.4 ± 0.8 

37.0 ± 0.4 

X 

18.9 ± 1.5 

-7.0 ± 1.5 

7.31-8.98 

-7.5 ± 0.7 ± 0.6 

-4.5 ±0.8 ± 0.3 

-8.5 ±0.8 ± 1.0 

16.5 ± 0.5 

Y 

38.1 ± 1.3 

16.2 ± 1.3 

11.31-14.84 

-15.1 ± 0.5 ± 2.0 

-3.3 ± 0.5 ± 0.2 

-0.4 ± 0.5 ± 0.2 

21.9 ± 0.4 

Z 

105.7 ± 2.0 

-85.7 ± 1.4 

17.58-23.06 

-11.4 ± 0.5 ± 

1.5 

-25.6 ±0.3 ± 2.6 

-18.1 ± 0.4 ± 3.2 

15.8 ± 0.5 


a: vsini = 34.0km /s l|Glebocki GnacinskilboOSh. b: Likely a Persei mem ber l|Hoogerwerf &: Blaauw||2OQ0h. c: Pleiades member _ 

llSarro et al.ll2Ql-3^ . usin^ = 11.0 zb 3.00 k m/s ijGlebocki &: Gnacinskill2005h. vsini = 1 2.2 di 1.0km/s iMermilliod. Mayor & Udr^l2009^ , 
d: vsini = 15.6 zb l.Okm/s llBiazzo et al.ll 20 l 4 ) . e: Parallax from ide Bruiine Sz Eilerj ||2012|) . 


Table 7. Properti es of the nearest known MGs within 100 pc. Column 2 is the number of MG objects with a measured parallax, from 
[Gagn^^^alJj 2^M ). The age rang es used for each MG are from th e following pub lications: T WA, Columba, Carina, Argus and ABDMG 
ijCagne et alJ2014 ) ; Octans-Near ijZuckerman et al.|[2013lb t? Cha l|Luhmanll2004l j and e Cha llMurphv . Laws on & Bess^lll20l'3^ . The age 
ranges for BPMG llBinks Sz Jeffriesll2014l : Malo et ^l2014al and Tucana-Horolo gium ('Tuc-Hor. iKraus et al.ll2014l) are solely from the 


‘li thium depletion bounda ry’ technique. UV WXYZ are from | Ga gne et ^ ||2014|) . except for Octans (calculated from the available data 
in IZuckerman et al.ll201^ and r] and e Cha llMurphv, Lawson &; Bessel]||2013h . Columns 7 and 8 correspond to the convergent points in 
right ascension and declination, respectively. 


Name 

N 

Age 

(Myr) 

Distance 

(pc) 

UVW 

cruvw 

(kms“^) 

XYZ 

O'XYZ 

(pc) 

acp 

D 

5cp 

(°) 

TWA 

12 

8-12 

42-92 

-10.53, -18.27, -5.00 
3.50, 1.17, 2.15 

12.17, -43.23, 21.90 
6.14, 7.30, 3.06 

180.87 

-79.93 

BPMG 

44 

21-26 

9-73 

-11.16, -16.19, -9.27 
2.06, 1.32, 1.35 

4.35, -5.82, -13.29 
31.43, 15.04, 7.56 

88.00 

-30.16 

ABDMG 

48 

70-120 

11-64 

-7.11, -27.21, -13.82 
1.39, 1.31, 2.26 

-2.25, 2.93, -15.42 
20.10, 18.97, 15.37 

92.89 

-47.73 

Tuc-Hor 

42 

39-43 

36-71 

-9.93, -20.72, -0.89 
1.55, 1.79, 1.41 

11.80, -20.79, -35.68 
18.57, 9.14,5.29 

116.28 

-28.80 

Columba 

20 

20-40 

35-81 

-12.24, -21.27, -5.56 
1.08, 1.22, 0.94 

-28.22, -29.74, -28.07 
13.68, 23.70, 16.09 

103.27 

-29.79 

Carina 

5 

20-40 

46-88 

-10.50, -22.36, -5.84 
0.99, 0.55, 0.14 

15.55, -58.53, -22.95 
5.66, 16.69, 2.74 

104.93 

-34.25 

Argus 

11 

30-50 

8-68 

-21.78, -12.08, -4.52 
1.32, 1.97, 0.50 

14.60, -24.67, -6.72 
18.60, 19.06, 11.43 

91.93 

-1.24 

Octans-Near 

14 

30-100 

24-98 

-13.04, -3.53, -11.24 
1.86, 2.17, 2.00 

8.07, 25.32, -48.97 
3.76, 11.81, 22.83 

60.12 

-4.74 

7] Cha 

4 

5-10 

~ 97 

-10.20, -20.70, -11.20 
0 . 20 , 0 . 10 , 0.10 

33.40, -81.00, -34.90 
0.40, 1.00, 0.40 

89.86 

-37.55 

e Cha 

35 

3-5 

100-120 

-10.90, -20.40, -9.90 
0.80, 1.30, 1.40 

54.00, -92.00, -26.00 
3.00, 6.00, 7.00 

92.48 

-35.13 
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Figure 7. Mk ns V — Ks CMD for the entire observed sample. 
Red lines indicate the range in for the likely-young sam¬ 

ple based on the assumed age range fo r each object. The model 
isochrones of ISiess. Dufour fc Forestinil ||2000|) are overplotted at 
5, 10, 30, 200 and 1000 Myr. The purple downwards-facing trian¬ 
gle is the absolute K magnitude for target P, measured from a 
trigonometric parallax. 


does not consider the possibility that our targets could be 
unresolved binaries which would enlarge our parallax errors. 

6.3 Comparison with MG space velocities 

To test if a star is kinematically matched in UVW to a 
known M G, a reduced fitti ng statistic is used, as de¬ 
scribed in I Shkolnik et all (|2012l '): 

-2 jU* -UMof (K - VMGf 

(7) 

With the exception of target P, the UVW midpoint 
values are used as inputs for UVW and half the line length 
is nsed for the errors. Table [8] provides the reduced Xt values 
for all of the likely-young sample. To interpret these results, 
given there are 3 degrees of freedom, a Xt value larger than 
3.78 rejects the null hypothesis at 99 per cent confidence. 
Any candidate/MG test where Xt > 3-78 results in rejection 
for membership of that MG. 

6.4 Predicted radial velocity comparisons with 
known MGs 

The Boettlinger diagrams require distances, which, with 
one exception for the likely-young sample, could only be 


10 

5 

^ I I I r 

I I ^ 

0 

-5 

M Octans-Near 

- ^ _ 

-10 

_ f Argus 


-15 

V C|ia ^ 

^TWA 
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-25 
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^ -In 
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Figure 8. Boettlinger diagrams in U — V and V — W space 
highlighting the space motions of the characteristically young, 
likely-single objects in this paper. Each object is connected by 
a solid purple line which indicates distance uncertainties due to 
the photometric parallax method described in mu Labels for the 
likely-young sample are the same as in Table [T] 


inferred using photometric parallaxes. Therefore a dis¬ 
tance independent criterion was used to strengthen the 
case for MG membership. We have recalculated the con¬ 
vergent p oints of each MG us i ng the membership lists pro¬ 
vided by iGagne e t, al.l llioidh. IZuckerman et al.l l 2013l (for 
Octans-Near) and Murphy, Lawson fc Besselll l2013r Tfor e 
and r] Cha). Given the UVW and convergent point of an 
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MG and the position of a candidate member, the predicted 
RV for a group member is then Vt cos A, where Vt is the 
total velocity of the group (\/U^ + V'^ + W^) and A is the 
angle between the convergent point and the target position. 

The average RV error for the likely-young sample is 
~ Ikrns”^ and the dispersion amongst MG members in 
UVW is < 3kms“^. Therefore true MG members are re¬ 
quired to have RVs that are within 5kms~^ of the predicted 
RV (ARV = |RV — VtcosA| < 5kms“^). Any candidates 
with ARV > 5kms“^ are rejected as MG members. Objects 
that have RVs within this threshold are flagged as possible 
members, although even if they satisfy this criterion, a dis¬ 
crepancy in distance, age or chemical abundances may result 
in a rejection. For this reason the technique can only be used 
to eliminate membership. The convergent p oint of each MG 
(following the prescription in Ijonesi Il97lli is given in Ta¬ 
ble [3 and values of ARV for each candidate/MG are given 
in Table |9l 

In 1)4.61 it was calculated that there is a 7.2 per cent 
chance that a star has 2 RV measurements within 5kms~^ 
but is still > 5kms“^ from the average centre-of-mass RV. 
This implies that out of a sample of 26 Li-rich stars likely 
to be single based on repeat RV measurements, ~ 2 may 
be shifted in UVW by ~ 5kms“^ which would potentially 
alter our assessment of their membership status for any MG. 

In principle we could perform a similar test using the 
tangential velocities, if the distance and proper motion are 
known accurately (i.e. JVtan < 5kms“^). In practice we 
can only do this on the one object - target P, which is within 
4.6kms“^ of ABDMG but > 5kms“^ for all other MGs. 


6.5 Combined MG membership criteria 

In this work a candidate only remains a possible member of 
a MG if it satisfies all of the following criteria: 

(i) The Xt value in Table[8]is no more than 3.78. 

(ii) ARV values in Table[9]are no more no than 5kms“^. 

(iii) The estimated age range for the candidates (provided 
in Table [5| must overlap with the age range of the MG, 
which are provided in a supplementary online table. 

In Table fTUl the outcome for each criteria is given. A i/ 
indicates a successful criterion match and an X represents a 
failure. 


6.6 Comparison with BANYAN 

iMalo et al.l ll2013ll and iGagne et al.l (l2014l l have presented 
tools for the Bayesian Analysis of Nearby Young Associ- 
atioNs (BANYAN I and BANYAN II respectively). These 
tools attempt to assign each star a probability of mem¬ 
bership for seven well-defined MGs (TWA, BPMG, Tuc- 
Hor, Columba, Carina, Argus, ABDMG) using sky posi¬ 
tion, proper motion, photometry and, where available, ra¬ 
dial velocities and parallaxes, as inputs. BANYAN II is an 
improvement over BANYAN I in terms of its treatment of 
uncertainties and a more realistic assessment of the likely 
field star contamin ation using the BesanQon Galactic model 
llRobin et aLllioi^ L 


We have used the provided web-based tool0 to test MG 
membership for our sample of 26 likely young objects. Using 
the BANYAN I technique, almost all our targets have zero 
or very low probability of membership for any of the tested 
MGs. In BANYAN II, the more realistic incorporation of 
the field population results in even lower probabilities; none 
of our young targets has a MG membership probability of 
more than 2 per cent. 

The results from BANYAN appear to contradict the 
assessments we have made in 46.51 The reason appears to 
be a philosophical difference in how MG membership is de¬ 
fined. BANYAN requires that new MG members have sim¬ 
ilar spatial coordinates to those members that are already 
known. The predominantly northern hemisphere locations 
of our targets mean that their XYZ positions are far from 
the centres of most of the MGs (as defined in BANYAN). 
Consequently, the membership probabilities are very low. 
However in this paper we choose to define members of MGs 
only on the basis of their youth and kinematics and so do 
not consider the BANYAN membership probabilities any 
further. 


7 INDIVIDUAL OBJECTS 

7.1 Individual objects linked to at least one MG 

Here we focus on objects that satisfied all 3 criteria in 46.51 
for at least one MG. The following section discusses any of 
these objects that have literature sources pertaining to MG 
membership and we compare them to our findings in this 
work. 

Target G - All three criteria are successfully matched with 
ABDMG. Although target G satisfies both the Xt s-nd 
RV criteria, we have not assigned it to BPMG because its 
age is measured between 30 and 150 Myr, too old to be a 
BPMG member (21-26 Myr, see Table [T]). The age of this 
object is more likely to be coeval with ABDMG (between 
70 and 125Myr). There is an uncertainty of ~ 4kms“^ in 
V velocity, which means the object could have kinematic 
matches to Tuc-Hor or TWA (but not BPMGL A litera¬ 
ture search for this object using the SIMBAEij database 
re veals that it was classified as a Pleiades me mber (Pels 22) 
by Ivan Leeuwen, Alphenaar fc Brandi (Il986l l and has since 
been used as a Pleiades member in at least 7 subsequent 
publications. Target G is very close to the centroid of the 
Pleiades cluster (a, 5 = 03h 47m 00s, -l-24d 07m 00s) and 
it has a RV consistent to within la, however the parallax 
range used in this work (8.78 — 11. 37 mas) does not encom¬ 
pass the Pleiades 17.34 ±0.06 mas. iMelis et al.ll20l3l . which 
casts doubt on its membership. If target G were an equal 
mass binary, an increase of 0.75 mag (half the flux) may re¬ 
solve the distance discrepancy. 

Target I - Although this K3 spectral-type object matches 
all three criteria for seven MGs (BPMG, Tuc-Hor, Columba, 
Carina, Octans, ry Cha and e Cha), there is a large uncer¬ 
tainty in V velocity. Figure [8] shows that membership with 
Tuc-Hor, Carina, Columba and Octans-Near cannot be ruled 


® available at http://www.astro.umontreal.ca/malo/baiiyaii. 
http://simbad.u-strasbg.fr/simbad/ 
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Table 8. Xx values for each candidate/MG match in the likely-young sample. 


Target 

TWA 

BPMG 

ABDMG 

Tuc-Hor 

Columba 

Carina 

Argus 

Octans-Near 

T) Cha 

e Cha 

A 

45.9 

22.2 

81.5 

38.3 

62.3 

186.8 

68.2 

3.2 

240.5 

45.7 

B 

63.1 

38.8 

116.9 

29.6 

73.4 

638.7 

38.8 

7.9 

11511.3 

64.6 

G 

16.1 

13.5 

18.2 

1.1 

6.1 

230.8 

27.9 

24.3 

1808.2 

13.0 

D 

62.1 

37.3 

111.8 

31.0 

67.2 

346.4 

28.5 

6.6 

1490.7 

61.9 

E 

10.8 

0.9 

13.0 

5.9 

3.2 

5.1 

17.4 

9.4 

13.4 

2.7 

F 

36.8 

17.3 

69.4 

22.6 

37.6 

92.4 

28.7 

3.7 

346.5 

33.1 

G 

10.4 

3.0 

1.4 

12.2 

13.6 

12.4 

55.8 

12.3 

65.6 

9.6 

H 

16.9 

4.4 

11.6 

5.5 

0.8 

3.2 

16.5 

12.8 

67.7 

6.5 

I 

14.1 

0.6 

3.0 

3.5 

1.3 

1.4 

19.1 

2.4 

1.6 

0.3 

J 

29.9 

21.8 

10.1 

9.5 

8.0 

12.1 

35.8 

39.9 

91.3 

13.7 

K 

50.9 

24.8 

96.3 

28.1 

55.5 

264.8 

20.9 

3.4 

1901.8 

47.0 

L 

18.3 

7.3 

10.1 

6.9 

2.2 

3.9 

11.3 

16.8 

176.9 

7.1 

M 

13.9 

3.4 

2.4 

10.8 

6.4 

8.4 

40.8 

16.4 

1.7 

1.3 

N 

265.1 

361.1 

146.2 

403.5 

602.5 

1599.4 

1039.7 

216.4 

5716.4 

342.3 

O 

16.8 

18.9 

11.6 

1.2 

8.9 

10.5 

45.6 

37.3 

69.7 

13.8 

P 

44.8 

29.2 

83.5 

24.1 

47.5 

132.2 

10.4 

10.8 

245.2 

47.4 

Q 

9.2 

5.9 

4.9 

5.0 

8.2 

5.6 

46.5 

16.6 

436.3 

11.7 

R 

59.8 

60.5 

77.1 

107.4 

144.6 

310.8 

206.9 

30.0 

347.7 

82.5 

S 

47.0 

23.9 

64.0 

56.8 

55.7 

103.5 

49.4 

27.2 

73.6 

50.4 

T 

25.5 

10.1 

39.1 

13.5 

12.9 

35.8 

3.2 

23.2 

161.6 

18.3 

U 

54.9 

29.2 

87.1 

38.6 

50.3 

90.2 

6.4 

9.6 

153.7 

51.9 

V 

28.0 

25.7 

11.2 

3.8 

9.9 

16.0 

43.0 

43.4 

55.8 

15.2 

W 

45.7 

33.5 

4.5 

29.4 

28.6 

74.6 

74.1 

46.2 

227.6 

16.9 

X 

52.4 

25.8 

96.3 

33.7 

65.7 

271.0 

41.8 

3.3 

832.5 

51.4 

Y 

68.8 

46.0 

124.9 

32.9 

81.3 

495.3 

27.3 

9.8 

3401.3 

73.2 

Z 

17.1 

5.8 

1.9 

9.2 

5.7 

5.6 

20.6 

15.6 

3.0 

3.0 


Table 9. ARV values for all MG/candidate comparisons 


Target 

TWA 

BPMG 

ABDMG 

Tuc-Hor 

Columba 

Carina 

Argus 

Octans-Near 

rj Cha 

e Cha 

A 

21.0 

5.5 

11.5 

13.8 

10.1 

11.5 

0.1 

4.3 

7.5 

7.7 

B 

17.2 

1.4 

6.4 

9.8 

5.8 

7.2 

3.6 

7.6 

3.0 

3.3 

C 

14.8 

4.5 

1.9 

2.4 

1.2 

0.5 

12.4 

12.8 

2.4 

2.5 

D 

23.6 

3.2 

9.0 

9.5 

5.9 

7.6 

5.1 

4.1 

5.0 

4.8 

E 

17.4 

1.0 

8.9 

2.9 

1.3 

3.4 

11.8 

10.3 

2.5 

1.8 

F 

22.1 

1.2 

6.6 

7.1 

3.5 

5.2 

7.4 

5.4 

2.8 

2.6 

G 

16.8 

4.5 

0.7 

1.1 

2.9 

1.2 

12.7 

9.6 

3.6 

3.7 

H 

23.2 

1.1 

6.0 

3.8 

1.1 

2.9 

9.3 

2.5 

2.4 

1.8 

1 

22.9 

0.9 

4.1 

3.1 

0.2 

1.8 

9.0 

0.9 

1.5 

1.0 

J 

19.2 

5.2 

9.8 

0.4 

1.4 

2.8 

5.5 

6.5 

6.6 

5.2 

K 

5.0 

7.2 

15.4 

0.9 

5.4 

6.1 

0.8 

4.7 

10.6 

9.2 

L 

10.9 

4.3 

2.7 

4.4 

0.4 

0.6 

7.2 

10.6 

4.1 

3.4 

M 

11.2 

2.8 

4.0 

9.2 

4.2 

4.1 

5.9 

4.7 

0.2 

0.9 

N 

59.9 

47.9 

40.6 

51.8 

47.2 

47.2 

48.5 

52.7 

44.6 

45.3 

O 

9.5 

9.1 

1.2 

2.4 

3.1 

2.4 

9.6 

21.3 

5.8 

5.5 

P 

1.6 

6.4 

10.5 

1.9 

3.1 

3.1 

3.5 

8.4 

8.3 

7.1 

Q 

14.4 

5.1 

5.6 

0.6 

4.1 

3.0 

10.2 

6.8 

5.9 

5.6 

R 

17.1 

4.5 

2.6 

8.3 

3.7 

3.8 

4.6 

9.2 

1.2 

1.9 

S 

26.8 

13.3 

6.4 

15.5 

11.4 

11.6 

12.2 

19.2 

10.1 

10.5 

T 

13.4 

3.9 

5.1 

5.6 

1.7 

1.4 

5.2 

11.3 

0.1 

0.5 

U 

8.9 

3.4 

11.0 

1.8 

5.7 

5.7 

4.1 

3.3 

7.0 

6.5 

V 

13.0 

5.0 

4.6 

5.3 

1.9 

1.6 

6.5 

13.9 

0.8 

1.3 

w 

17.5 

13.6 

3.0 

11.5 

9.4 

8.7 

16.1 

24.7 

9.1 

9.5 

X 

2.5 

8.3 

9.1 

14.4 

13.7 

13.0 

15.5 

1.7 

9.1 

10.1 

Y 

9.7 

5.6 

12.1 

14.1 

11.9 

12.7 

5.0 

5.9 

8.2 

8.7 

Z 

2.7 

1.7 

4.0 

7.1 

4.5 

5.2 

2.0 

12.8 

0.6 

1.1 


out, however a more accurate parallax is required to iden¬ 
tify which group it may belong to. The object has an Li 
EW (> 300 m A), large enough to set an age < 30Myr. 
iLi Sz HtJ lll998l i identify this object as a weak-lined T-Tauri 
star (WTTS) in the surrounding area of the Taurus-Auriga 
region. It is classified as a T-Tauri star in the pre-main se¬ 


quence catalog developed bv iDucourant et ahl ll2005h . how¬ 
ever no excess emission indicative of circumstellar material 
is observed at near-infrare d wavelength s. The distance to 
Taurus (140pc, ~ 7.1 mas, [Chd^[2003) is within our pho- 
tom etric parallax range of 6.46 — 12.74 mas. 

ISantos et al.1 ll2008ll located the object in the Taurus 
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Table 10. The full list of candidate/MG criteria matches for the likely-young sample. A ^ marks a successful match and a 1C denotes a 
failed criteria. The columns labelled T, 2 and 3’ are the criteria listed in 163] An object was only considered a potential member of an 
MG if all three were satisfied. 


TWA 
Target 1 2 

3 

BPMG 

12 3 

ABDMG 
12 3 

Tuc-Hor 

12 3 

Columba 

12 3 

Carina 

1 2 

3 

Argus 
1 2 

3 

Octans-Near 

12 3 

r) Cha 
1 2 

3 

e Cha 

1 2 

3 

A 

K 

X 

X 

X 

X 

X 

X 

X 

✓ 

X 

X 

✓ 

X 

X 

✓ 

X 

X 

✓ 

X 

✓ 

✓ 

✓ 

✓ 

✓ 

X 

X 

X 

X 

X 

X 

B 

K 

X 

X 

X 

✓ 

X 

X 

X 

✓ 

X 

X 

✓ 

X 

X 

✓ 

X 

X 

✓ 

X 

✓ 

✓ 

X 

X 

✓ 

X 

✓ 

X 

X 

✓ 

X 

C 

X 

X 

X 

X 

✓ 

X 

X 

✓ 

✓ 

✓ 

✓ 

✓ 

X 

✓ 

✓ 

X 

✓ 

✓ 

X 

X 

✓ 

X 

X 

✓ 

X 

✓ 

X 

X 

✓ 

X 

D 

X 

X 

X 

X 

✓ 

X 

X 

X 

✓ 

X 

X 

✓ 

X 

X 

✓ 

X 

X 

✓ 

X 

X 

✓ 

X 

✓ 

✓ 

X 

✓ 

X 

X 

✓ 

X 

E 

X 

X 

X 

✓ 

✓ 

X 

X 

X 

✓ 

X 

✓ 

✓ 

✓ 

✓ 

✓ 

X 

✓ 

✓ 

X 

X 

✓ 

X 

X 

✓ 

X 

✓ 

X 

✓ 

✓ 

X 

F 

X 

X 

X 

X 

✓ 

X 

X 

X 

✓ 

X 

X 

✓ 

X 

✓ 

✓ 

X 

X 

✓ 

X 

X 

✓ 

✓ 

X 

✓ 

X 

✓ 

X 

X 

✓ 

X 

G 

X 

X 

X 

✓ 

✓ 

X 

✓ 

✓ 

✓ 

X 

✓ 

✓ 

X 

✓ 

✓ 

X 

✓ 

✓ 

X 

X 

✓ 

X 

X 

✓ 

X 

✓ 

X 

X 

✓ 

X 

H 

X 

X 

X 

X 

✓ 

X 

X 

X 

✓ 

X 

✓ 

✓ 

✓ 

✓ 

✓ 

✓ 

✓ 

✓ 

X 

X 

✓ 

X 

✓ 

✓ 

X 

✓ 

X 

X 

✓ 

X 

I 

X 

X 

✓ 

✓ 

✓ 

✓ 

✓ 

✓ 

X 

✓ 

✓ 

✓ 

✓ 

✓ 

✓ 

✓ 

✓ 

✓ 

X 

X 

✓ 

✓ 

✓ 

✓ 

✓ 

✓ 

✓ 

✓ 

✓ 

✓ 

J 

X 

X 

X 

X 

X 

X 

X 

X 

✓ 

X 

✓ 

✓ 

X 

✓ 

✓ 

X 

✓ 

✓ 

X 

X 

✓ 

X 

X 

✓ 

X 

X 

X 

X 

X 

X 

K 

X 

✓ 

X 

X 

X 

X 

X 

X 

✓ 

X 

✓ 

✓ 

X 

X 

✓ 

X 

X 

✓ 

X 

✓ 

✓ 

✓ 

✓ 

✓ 

X 

X 

X 

X 

X 

X 

L 

X 

X 

X 

X 

✓ 

X 

X 

✓ 

✓ 

X 

✓ 

✓ 

✓ 

✓ 

✓ 

X 

✓ 

✓ 

X 

X 

✓ 

X 

X 

✓ 

X 

✓ 

X 

X 

✓ 

X 

M 

X 

X 

X 

✓ 

✓ 

X 

✓ 

✓ 

✓ 

X 

X 

✓ 

X 

✓ 

✓ 

X 

✓ 

✓ 

X 

X 

✓ 

X 

✓ 

✓ 

✓ 

✓ 

X 

✓ 

✓ 

X 

N 

X 

X 

X 

X 

X 

X 

X 

X 

✓ 

X 

X 

✓ 

X 

X 

✓ 

X 

X 

✓ 

X 

X 

✓ 

X 

X 

✓ 

X 

X 

X 

X 

X 

X 

O 

X 

X 

X 

X 

X 

X 

X 

✓ 

✓ 

✓ 

✓ 

✓ 

X 

✓ 

✓ 

X 

✓ 

✓ 

X 

X 

✓ 

X 

X 

✓ 

X 

X 

X 

X 

X 

X 

P 

X 

✓ 

X 

X 

X 

X 

X 

X 

✓ 

X 

✓ 

✓ 

X 

✓ 

✓ 

X 

✓ 

✓ 

X 

✓ 

✓ 

X 

X 

✓ 

X 

X 

X 

X 

X 

X 

Q 

X 

X 

X 

X 

X 

X 

X 

X 

✓ 

X 

✓ 

✓ 

X 

✓ 

✓ 

X 

✓ 

✓ 

X 

X 

✓ 

X 

X 

✓ 

X 

X 

X 

X 

X 

X 

R 

X 

X 

X 

X 

✓ 

X 

X 

✓ 

✓ 

X 

X 

✓ 

X 

✓ 

✓ 

X 

✓ 

✓ 

X 

✓ 

✓ 

X 

X 

✓ 

X 

✓ 

X 

X 

✓ 

X 

S 

X 

X 

X 

X 

X 

X 

X 

X 

✓ 

X 

X 

✓ 

X 

X 

✓ 

X 

X 

✓ 

X 

X 

✓ 

X 

X 

✓ 

X 

X 

X 

X 

X 

X 

T 

X 

X 

X 

X 

✓ 

X 

X 

X 

✓ 

X 

X 

✓ 

X 

✓ 

✓ 

X 

✓ 

✓ 

✓ 

X 

✓ 

X 

X 

✓ 

X 

✓ 

X 

X 

✓ 

X 

U 

X 

X 

X 

X 

✓ 

X 

X 

X 

✓ 

X 

✓ 

✓ 

X 

X 

✓ 

X 

X 

✓ 

X 

✓ 

✓ 

X 

✓ 

✓ 

X 

X 

X 

X 

X 

X 

V 

X 

X 

X 

X 

✓ 

X 

X 

✓ 

✓ 

✓ 

X 

✓ 

X 

✓ 

✓ 

X 

✓ 

✓ 

X 

X 

✓ 

X 

X 

✓ 

X 

✓ 

X 

X 

✓ 

X 

w 

X 

X 

X 

X 

X 

X 

X 

✓ 

✓ 

X 

X 

✓ 

X 

X 

✓ 

X 

X 

✓ 

X 

X 

✓ 

X 

X 

✓ 

X 

X 

X 

X 

X 

X 

X 

X 

✓ 

X 

X 

X 

X 

X 

X 

✓ 

X 

X 

✓ 

X 

X 

✓ 

X 

X 

✓ 

X 

X 

✓ 

✓ 

✓ 

✓ 

X 

X 

X 

X 

X 

X 

Y 

X 

X 

X 

X 

X 

X 

X 

X 

✓ 

X 

X 

✓ 

X 

X 

✓ 

X 

X 

✓ 

X 

✓ 

✓ 

X 

X 

✓ 

X 

X 

X 

X 

X 

X 

Z 

X 

✓ 

X 

X 

✓ 

X 

✓ 

✓ 

✓ 

X 

X 

✓ 

X 

✓ 

✓ 

X 

X 

✓ 

X 

✓ 

✓ 

X 

X 

✓ 

✓ 

✓ 

X 

✓ 

✓ 

X 


SFR and derive chemical abundances. To test for abundance 
matches with MGs, the mean and standard deviations of 
[Fe/H], [Si/ H] and [Ni/H] for mem bers of the MGs in tables 
2 and 3 in IViana Almeida et al.l (l2009l i were compared to 
target I. The abundances are consistent with any of the con - 
sidered MGs to within 1 a uncertainties. iKazzo et al.l(l2012l l 
identify this target as part of a 2 — 10 Myr stellar aggregrate 
within 15° of Taurus-Auriga, which they name ‘X-Clump 


0534-b22’. 

Target Z - This target is matched with ABDMG and 
membership is rejected for TWA because target Z is too 
old to be coeval with TWA. The Li EW/F — Ks of the 
target sets an a ge ra nge of 30 — 150 Myr (see Figure [S]). 
iLi, Hu CherJ 1 2000fl ident i fy th e object as a possible 
WTTS and Ducourant et al.l (l2005l i include it as a member 
in the PMS stars proper-motion catalog. The UVW values 
all lie outside the ABDMG box in Figure [Hi however, the 
UVW uncertainties from parallax estimation are large. A 
more precise age and parallax are necessary to test if it is 
an ABDMG member. 


7.2 Individual objects not linked to any MG 

Target B - At a spectral-type of G1 and Li EW of 154 mA 
this target is assigned an age range of 30 — 200 Myr, sup¬ 
ported by a lower age limit of 20 Myr using Ha measure¬ 
ments. In terms of UVW velocity, it appears consistent in 
U — V with Octans-Near, however the W v elocity of the 
target is ~ lOkms”^ larger than Octans-Near. Ihhvata et aH 
(Il995ll measure a v sin i value of 34 km s ^ (no error bars 


provided), which is 1.5 kms ^(3 (t) slower than the measure¬ 
ment in this work. 

Target D - The age measured for the target is 30 — 200 Myr 
based on Li measurements. This age is supported by a lower 
age limit of 20 Myr set by Ha measurements. It is in the 
range of Octans-Near members in terms of U and V velocity, 
but is ~ 10kms“^ larger in W. It is flagged as a c andidate 
WTTS in Taurus-Auriga region in lLi fc M lll998ll and our 
photometric parallax range is consistent with this idea. 

Target F - Age estimates are 30 — 200 Myr from Li measure¬ 
ments, and > 25 Myr based on Ha measurements. It appears 
consistent with the U and V velocities of Octans-Near, but 
is around 5kms“^ larger in W. 

Target J - An age range of 30 — 150 Myr was estimated 
from Li measurements, consistent with a lower age limit of 
20 Myr set by Ha. It is > 5kms“^ in U, V and W from any 
of the 10 MGs considered. 

Target N - Although target N has an Li EW of 86 mA, 
because it is a late-K dwarf it is estimated to have an age 
of 30 — 150 Myr, with an upper limit of 189 Myr set by Ha 
measurements. The W velocity is matched to a few MGs, 
however, both the U and V are at least ~ 15kms“^ from 
any MG. 

Target P - Using Li, an age of 30 — 150 Myr was estimated 
and Ha measurements suggest an age of 65 — 125 Myr. It 
has a U,V and W v elocity at least lOhms”*^ from any MG. 
IXing fc Xind ll2012h include this target as a young K4V ob¬ 
ject (we calculate K5) with a Li EW of 17 mA which is 
almost 150 mA lower than the measurem ent in thi s work . 
It is listed as a possible visual binary in iMakarp-^ ll2003ll , 
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however, two consistent RV measurements in this analysis 
suggest the object is uot a dose biuary. 

Target Q - The age range of 30 — 150 Myr may be more 
consistent with the pattern observed in IC 2602 given the 
spectral-type K5 and Li EW of 174 mA. The Ha age is es¬ 
timated as 40-65 Myr. The U and V velocities are consis¬ 
tent with ABDMG, however its W velocity is too large by 
~ inkms~A 

Target R - An age of 30 — 150 Myr constrained by Li mea¬ 
surements is supported by an age of > 20 Myr from Ha. It 
is at least lOkms”^ from any MG in either U, V or W. 
Target S - Because target S has a late-K spectral-type, its 
Li EW (135 mA) leads to an age of 30 — 150 Myr, supported 
by a Ha age of < 230 Myr and an R sin i age of < 80 Myr. 
It is > 5kms“^ from any MG in U, V or W. 

Target T - Its Li EW is only 62 mA, however, as it is a late- 
K type, it was estimated to have an age of 100 — 200 Myr. 
This is too old to be a BPMG member. An limit of < 
120 Myr is provided because Ha is in emission. There are 
no matches < 5kms“^ in all UVW velocities for any MG. 
Target U - This target has an age range of 30 — 150 Myr, 
however given the large Li EW (278 mA), this is more likely 
towards the younger end of this range. It has a Ha age older 
than 35 Myr. There are no matches within 5 km s“^ for U, V 
or W for any MG and BPMG membership can be ruled out 
based on both the kinematics and age. 

Target V - The age range is estimated as 30 — 150 Myr 
from Li and < 50 Myr based on Ha emission. Target V is 
not within 5kms~^ in U, V or W for any MG. 

Target W - This object satisfies the RV criteria in 46.41 
for ABDMG (see Table [O]) and its location in U — V and 
V — W velocity space suggests that it is similar to AB¬ 
DMG in V and W but is more than 5kms“^ from AB¬ 
DMG in U velocity and therefore fails membership test 1, 
although only marginally (xt = 4.5, the membership cri¬ 
teria is Yf < 3.78). The U, V and W are precise to within 
Ikrns”^ in each velocity coordinate. The Li EW measure¬ 
ment of 222 mA at an early-K spectral-type results in an 
age range of 30 to 150 Myr and although there is signifi¬ 
cant scatter amongst IC 2602 and Pleiades objects at this 
spectral-type, this target is more consistent with the lower 
envelope of IC 2602. The Ha age range is 35 — 55 Myr and 
an age of < 80 Myr is im plied by Rsini. 

iFrasca et al.l (l201ll ~l identify differential rotation based 
on data from the Kepler mission. Their Li EW measure¬ 
ment is ~ 50 mA larger and their I.OA Ha EW is certainly 
not consistent with the 0.1 A absorption line measured in this 
work. The v sin i measurements match to Icr and their calcu¬ 
lated inclination angle, i (~ 70°), provides a stellar radius of 
O.937?0, which from Figure [6] suggests an age of25 — 80 Myr, 
where the sin i ambiguity has been resolved. This is at least 
40 Myr younger than ABDMG (if indeed ABDMG is coeval 
with the Pleiades), making membership for ABDMG un¬ 
likely. The authors also note that the object is flagged as an 
equal-magnitude vis ual binary system in the WPS Catalog 
(WDS 19251-k4431, IWorlev fc Douglass! 1 1997l'l a nd the Ty¬ 


cho Double Star Catalog llFabricius_et^_aL 2002) . Base d on 
the orbital properties calculated by Frasca et al. (l201ll) . if it 


were an equal-mass binary it would be expected to have RV 
variations > 5kms ~^ over 3 — 4 years. The RVs reported by 
iFrasca et al.l (l201ll) were from 3 separate observations be¬ 
tween 2007 and 2009 (—33.1,—33.1,—32.2kms“^) and the 


measurement of —33.5 ±0.6kms“^ in this work carried out 
in June 2011 suggests that the object is not in a equal-mass 
visual binary system. 

Target Y - This target has a spectral-type K4 and a Li 
EW of 141 mA, from which an age range of 30 — 150 Myr 
was estimated, and a Ha emission line implies an age of 
< 100 Myr. It is within the range of Octans-Near in U and 
V velocity, but has a W velocity ~ 5kms“^ larger than 
Octans-Near. 


7.3 Connection with the Octans-Near association? 

In a Hipparcos-based survey, IZuckerman et all (l2013l) re¬ 
port 14 star systems with spectral-types between G5 and 
AO and distances ranging from 24 to 98 pc that have ages 
between ~ 30 — 100 Myr and Galactic space velocities sim¬ 
ilar to the mean space motion of Octans (l/VlFoctans = 
-14.5 ± 0.9, -3.6 ± 1.6, -11.2 ± 1.4kms~^ l/ VlFsampie = 
-13.0 ± 1.9, -3.5 ± 2.2, -11.2 ± 2.0kms”^). [Torres et all 
ll2008l) measure the distance of Octans to be 141 ± 34 pc and 
given that the 14 star systems in IZuckerman et al.l (l2013l ) are 
closer than the Octans association (all within 100pc), the 
authors describe this stellar aggregate as the ’Octans-Near’ 
group. 

In 46.51 four objects - A, I, K and X - pass tests for 
Octans-Near membership, however, target I is unlikely to be 
an Octans-Near member because the UVW error bars are 
so large that the fact it agrees with Octans-Near is not re¬ 
markable (see ani. Whether these objects are connected to 
either the Octans, the Octans-Near or not, there neverthe¬ 
less still exists a sub-grouping of 7 objects in the sample (A, 
B, D, F, K, X and Y) which are possibly coeval. All 7 ob¬ 
jects could be consistent with a common age of between 30- 
150 Myr, which is also consistent with Octans-Near. These 
have UVW values of 1/ = —10.5 ± 2.7(±1.8); V = —4.3 ± 
1.0(±1.4); VF = —4.9 ± 3.6(±1.3) kms“^, which is close to 
Octans/Octans-Near but not to any other known MG (error 
bars are the standard deviations and the values in paren¬ 
theses are the average UVW uncertainties). Of these seven, 
none are in similar in terms of XYZ positions, although five 
of them, A, B, D, F and Y, have sky positions within 3 hours 
of right ascension of one another and 10 degrees in declina¬ 
tion. These are in the vicinity of the Pisces constellation, 
therefore if these did constitute a MG and were not associ¬ 
ated with Octans-Near, we propose naming the ensemble as 
“the Pisces moving group”. 

In principle, if the 7 kinematically consistent objects 
identified here were born from the same molecular cloud, 
then given their UVW XY Z coordinates, it would be possi¬ 
ble to trace-back their motions under a reasonable gravita¬ 
tional potential and identify a time at which all these objects 
occupied a minimum. However, since our objects all have 
designated ages older than 30 Myr and lack trigonometric 
parallaxes, a kinematic traceback analysis is premature: it 
would only take 5 Myr for the current velocity uncertainties 
to produce a spatial dispersion equal to the current standard 
deviation of stellar positions. 
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8 CONCLUSIONS 

From an initial sample of 146 spectroscopically observed, 
FGK stars with short rotation periods from the Super WASP 
All Sky Survey, 26 were found to be Li-rich and there¬ 
fore probably young by comparison with the photospheric 
Li abundance patterns seen in open clusters. Duplicate RV 
measurements suggest these are unlikely to be close bina¬ 
ries. Five objects were G-stars and the remaining 21 have 
spectral types from KO to K6. 

Twelve objects from the initial observed target sample 
had 2 or more RV measurements that varied by > 5kms~^ 
on the timescale of the observing run and are almost cer¬ 
tain to be close, tidally-locked binary systems. This fraction 
of close binary systems is much higher than expected from 
a random sample of field stars, presumably because of our 
initial selection based on rapid rotation. If a signihcant Li 
line was not observed at the telescope, then usually a repeat 
measurement was not made. Of the remaining Li-depleted 
objects, for which only a single RV measurement was ob¬ 
tained, we anticipate that a large fraction will turn out to 
be close binary systems. 

Typical ages of members in the likely-young sample are 
30-200 Myr and have distances of 30-180 pc. Ages are pri¬ 
marily determined from lithium, but ages derived from rota¬ 
tion, activity and projected radius are usually in agreement, 
but might be compromised by binarity. Lower-mass, nearby, 
young stars are favourable as targets to follow-up with high- 
resolution imaging to detect the presence of brown-dwarf bi¬ 
naries, or even exoplanets. There are 16 objects in the likely- 
young sample with ages < 200 Myr, later than KO and within 
100 pc that would be ideal candidates for imaging surveys. 

Fifty-four per cent of the observed sample had rotation 
periods < 2 days (77 out of 146) and there were fourteen 
objects in the likely-young sample that had rotation peri¬ 
ods less than 2 days. The sample may not be complete since 
sometimes even young stars do not exhibit significant rota¬ 
tional modulation. In terms of the distribution of rotation 
periods, the likely-young sample and parent sample have 
similar medians, therefore one would not gain by being more 
restrictive in choosing a rotation period cut-off. 

Eleven targets are identified which satisfy kinematic 
and age criteria with at least one MG and fifteen fail at least 
one criteria for membership of any of the MGs considered in 
this analysis. There is a low probability of the RVs being sig¬ 
nificantly in error due to binary motion (see 44.61) . therefore 
it is unlikely that many of these classifications would change 
with more extensive RV observations. There is tentative ev¬ 
idence for a sub-grouping of seven objects that appear to be 
a northern hemisphere counterpar t to, but not identical to , 
the ‘Octans-Near’ MG proposed bv IZuckerman et al.N2013li . 
Regardless of whether some of these objects are connected 
to Octans-Near or not, a grouping of stars exist around 
U = -10.5 ± 2.7; V = -4.3 ± 1.0; VF = -4.9 ± 3.6kms"i 
which we tentatively label as “the Pisces moving group”. 

In light of results from the analysis of these 146 ob¬ 
jects, additional substantial surveys of MG candidates in the 
northern hemisphere are desirable. The kinematically unbi¬ 
ased search mechanism was ~ 18 per cent efficient at detect¬ 
ing likely-single G and K-type stars younger than 200 Myr 
and can be used to identify stars that are not associated 
with previously studied MGs. More widespread searches for 


MG candidates in the northern hemisphere may provide ad¬ 
ditional useful samples of young stars, and could unveil im¬ 
portant kinematic substructure in the Solar neighbourhood. 
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